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PART  I 


AH  EXrERi MEN  < AL 


3 j UDY  Or  rArcFiAL  AND  COMrLtic  HLM  builinl 


WITH  FORCED  CONVECTION 


SUMMARY 


The  conventional  laboratory  equipment  for  studying  boiling  beet  transfer  was  modified  no 
that  it  could  be  operated  in  a stable  manner  with  subcode'"  liquids  in  the  regions  cf  .,  .'leate,  parti"! 
film,  and  complete  film  boiling.  The  apparatus  employed  a econdary  stabilizing  fluid  which  flowed 
through  the  inside  of  an  electrically-heated  stainless-steel  tube  while  the  tent  imid  w.-a.  flov.1.7 
through  an  annulus  formed  by  the  tube  and  a Pyrex  jacket.  The  stabilizing  fluid  absorbed  the  excess 
heat  which  could  not  be  transferred  to  the  test  fluid.  This  arrangement  allowed  the  apparatus  to 
operate  safely  in  all  three  boiling  regions  since  the  total  heat  transferred  to  the  stabilizing  fluid 
and  tbe  test  fluid  was  a mcnotonically  increasing  function  of  wall  temperature  up  to  tbc  melting 
point  of  tbe  wall. 


T’l.n  am.! 

A 44V. 


mp’icny  of  eiee  ni  healing,  w*»s  used  to  study  boiling 
in  distilled  water  which  was  flowing  at  various  dociti^-,  pressures,  and  temperatures  in  an  annulus. 
The  results  of  this  investigation  and  a description  of  tbe  appuratus  are  presented  together  with  an 
approximate  method  for  calculating  the  heat  transfer  in  the  complete  film  boiling  region. 


High-speed  motion  pictures  that  were  taken  of  the  degassed  water  boiiing  on  the  electrieally- 
hsatsd  tube  showed  the  types  of  vapor  formation  in  the  three  boiling  regions.  The  mechanism  cf  tran- 
sition from  uucleete  to  partial  film  and  finally  to  complete  film  boil'”'  *s  the  wall  temperature  was 
increased  is  discussed. 


I. 


INTRODUCTION 


The  heat  transfer  from  a hot  metal  wall  to  a boiling  liquid  has  been  the  subject  of  extensive 
studies  which  started  with  the  work  of  Mosciki  and  Broder  (Cf.  R'  1).  Little  progress  was  made 
until  the  need  for  high  rates  of  heat  transfer  in  nuclear  reactors  and  rocket  motors  added  to  the 
desirability  of  obtaining  a fuller  understanding  of  the  boiling  heat  transfer  mechanics.  These  recent 
incentives  have  reunited  in  contributions  to  the  field  of  nucleate  boiling  by  many  iuveatLj'ors  in- 
cluding! McAdams  (Cf.  Ref.  2),  Trarnontini  (Cf.  Ref.  3),  and  Gunther  (Cf.  Ref.  4). 

It  was  noi  possible  for  these  investigators  to  make  a o.-  led  study  of  the  boiling  mechanism 
in  all  regions  because  of  the  .-i-  bility  of  er'sting  equipment  t-  -naii.1  stable  operation  near  the 
peak  heat  flux  and  partial  film  boiling  regions.  The  following  discussion  of  the  three  boiling  regions 
brings  out  this  inherent  difficulty. 

Figure  1 shows  a typical  heat  transfer  curve  in  which,  for  simplicity  of  the  discussion,  heat 
flux  is  plotted  against  wail  temperature  instead  of  the  customary  parameter  of  temperature  difference 
between  the  hat  wall  and  the  fluid.  The  shape  of  the  curve  is  unaffected  by  this  substitution  for 
constant  fluid  temneratnre.  nressure.  and  veloeitv. 

The  region  A to  B in  Figure  i represents  the  nonboiling  line  for  either  forced  or  free  convection. 
When  the  wall  temperature  increases  above  the  liquid  saturation  temperature,  such  as  at  point  B,  the 
liquid  against  the  hot  well  boils  at  small  nuclei.  If  the  bulk  of  the  liquid  is  below  saturation  temperature, 
the  bubbles  grow  to  a maximum  radius  and  then  collapse.  Since  the  boiling  occurs  from  small  active 
nuclei  at  the  surface  of  the  hot  wall  in  the  region  i’rom  B to  C,  this  type  of  boiling  is  called  surface  or 
nucleate  hailing.  High-speed  motion  pictures  of  nucleate  boiling  were  taken  with  the  fields  o?  view 
shown  in  Figure  2,  Several  frames  of  nucleate  boiling  are  shown  in  Figure  3. 
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/ s the  wall  temperature  is  increased  to  point  C in  Figure  1,  the  population  of  nucleate 
bubbles  increases  until  thev  become  so  crowded  that  teveiul  of  the  individual  bubbles  coalesce 
to  form  a local  blanket  of  vapor.  These  vapor  blankets  cover  localised  areas  of  tbr  hot  wall  for 
periods  of  about  ten  to  twenty  timers  the  lifetime  of  a nucleate  boiling  bubble  befere  condensing 
is  the  main  streK.m. 

As  the  wall  temperature  is  inorea.'ieri  still  further,  such  as  from  C!  to  D of  h igure  1,  the 
wall  area  covered  by  the  local  fil.n  ult  aVeis  increases  and  causes  the  over-all  heat  flux  to  de- 
crease. The  region  between  C and  D is  called  Partial  film  boiling,  because  of  the  intermittent 
coverage  of  the  wai".  by  the  vapor  film.  Figure  4 shows  one  cycle  of  partial  film  boiling. 

At  point  D of  Figure  1 the  local  vapor  blankets  coalesce  to  form  a continuous  vapor  film 


This  is  the  complete  film  boiling  region.  In  this  region  t 


the 


to  the  liquid  principally  by  conduction  across  the  laminar  vapor  film.  A further  increase  of  wall 
temperature  causes  a small  increase  in  vapor  film  thickness.  The  thickness,  however,  increases 
at  a lower  rate  the.n  the  temperature  difference  ATsa,,  so  that  the  heat  flax  increases.  At  higher 
wall  temperatures  the  film  gets  so  thick  that  it  is  no  longer  laminar.  The  roughness  of  the  film, 
coupled  with  radiation,  makes  the  heat  transfer  increase  more  rapidly  with  an  increase  in  *”*11 
temoerature.  This  increase  results  in  the  concave  upward  curve  E to  F in  Figure  1.  The  rough 
i il i:t  i»  s It >.■  „ n in  vie »v  u of  Figure  b.  tough  film  would  occur  eitucr  Fit  fiigh  wall  temperature 

or  on  a very  long  heating  section.)  The  vapor  film  is  also  visible  it  view  b of  Figure  5 which 
shows  the  film  as  it  is  being  detached  from  the  heating  tube  by  the  sharp  corner  of  the  copper  bus 
bar. 

In  addition  to  the  three  regions  of  boiling  heat  transfer,  the  boiling  process  may  be  divided 
into  two  general  categories:  one,  where  the  bulk  temperature  or  the  li-quid  is  below  its  saturation 
temperature  (hest  transfer  to  subcooled  liquids),  and  the  other,  where  the  liquid  bulk  temperature  is 
at  its  saturation  temperature  ( bulk  boiling).  The  fluids  used  in  the  regenerative  cooling  of  a rocket 

. » _ i - r i •!?_  • i _ . i _ J l : _ • J — l .l*  i.' ? - L . : I 

lliUkur  arc  ail  example  ui  ooiuttg  hi  a i:\jiiiti.  niicican  uic  gcuciauuu  ui  diccmii  tu  a uuu.i 

is  an  example  *r  bulk  boiling. 

It  is  now  possible  to  discuss  the  difficulty  of  measuring  the  heat  flux  in  the  three  boiling 
regions.  The  conventional  laboratory  apparatus  employed  for  the  study  of  heat  transfer  has  a high 
resistance  wall  (usually  a thin  metal  strip  or  tube)  which  is  heated  electrically.  Reference  to  the 
solid  curve  in  Figure  6,  which  represents  heat  transfer  to  a subcooied  iiquid,  shows  that  up  to 
point  A any  increase  in  electrical  power  supplied  to  ihc  melts'  wall  may  be  absorbed  uv  i’ne  liquid. 
Thus  the  convection  and  nucleate  boiling  regions  may  be  studied  without  difficulty.  However,  if 
the  power  is  increased  above  point  A.  the  wull  temperature  must  jump  to  poiut  R before  the  additional 
^iwer  may  be  transferred  to  the  liquid.  At  point  13,  the  heat  transfer  surface  temperature  required  for 
Tiost  liquids  is  above  the  melting  point  of  the  metallic  wall,  and  failure  occurs. 

It  has  been  possible  to  study  the  nucleate-boiling  and  complete-film-boiling  regions  for  a 
hulk  boiling  liquid  as  is  represented  by  the  dashed  line  of  Figure  6.  The  peak  heat  flux  shown  a* 
point  A is  considerably  lower  than  for  the  same  iiquid  when  subcooled.  It  is  seen  to  be  possible 
also  to  transfer  the  peak  quantity  of  heat  in  the  complete-film-boiling  region  without  exceeding  the 
melting  point  of  the  wall. 

Figure  7 shows  Nukiyania’s  curve  (Cf.  Ref.  5)  of  the  heat  transferred  to  water  from  a 
G.GGC5-inch  platinum  wire  in  pool  bulk  boiling.  T’nr  peak  heat  transfer  shown  at  point  A is  suffi- 
ciently Sow  to  allow  the  wall  teraperatuie  to  remain  below  the  melting  point  in  the  film-boiling 
region  and  still  transfer  more  heat  than  at  the  peak  point.  It  is  thus  possible  to  study  t"h*  nucleate- 
boiling  and  complete-film- boiling  regions  tor  a stagnant  pooi  if  the  water  is  at  .he  saturation  tem- 
perature. Nukiyama’s  apparatus  could  not  be  used  to  investigate  the  partial-film-boiling  region. 

For  the  case  where  the  water  is  not  at  the  saturation  temperature,  there  are  (our  possible 
solutions  that  would  still  allow  the  study  to  be  extended  into  the  compie'.e-fiim-boiling  region. 
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One  method  would  be  to  control  the  wall  terr.pcrnture  rather  than  the  power  supplied  to  the  wall.  In 
this  case,  instead  of  using  electrical  heating,  tonUoi  could  be  accomplished  by  using  condensing 
vwnors  nr  fnsi-finwing  hot  liquids  that  would  supply  the  heat  to  the  wall.  This  arrangement  is  not 
desirable  since  it  has  many  serious  difficulties  in  control,  installation,  and  measurement  of  the 
amount  of  heat  supplied  to  the  wall,  A second  method  for  eliminating  burnout  but  which  still  retains 
ilic  desirable  electrical  Heating  would  be  to  make  the  wall  of  n high  melting  point  material.  However, 
lor  any  eiectricai  conductor,  this  method  wou  d prevent  burnout  only  over  a narrow  range  cf  liquid 
properties.  A third  method  of  preventing  burnout  would  be  to  decrease  the  electrical  powf  .wpe’ied 
to  the  wall  immediately  after  it  is  ruised  above  point  A of  Figure  6.  If  the  power  were  carefully  con- 
trolled, it  would  be  possible  to  jump  from  point  A to  some  other  point  in  the  film  boiling  region  that 
required  a wall  temperature  lower  than  point  B.  This  method  requires  elaborate  electtonic  equipment 
that  also  would  not  be.  satisfactory  except  for  limited  ranges  of  liquid  properties.  Of  these  three 
methods,  only  the  first  one  could  operate  in  the  partial-film-boiling  regiou. 

A fourth  method  will  he  considered  which  makes  it  possible  to  retain  the  convenience  of 
electrical  heating  and  still  maintain  stable  operation  in  ail  three  boiling  regions  for  a wide  range  of 
liquid  variables.  As  has  been  seen,  the  difficulty  with  the  conventional  equipment  which  uses  elec- 
trical heuting  is  that  a small  increase  in  power  above  the  peak  point  causes  the  wail  temperature  to 
rise  above  the  melting  point,  if  the  heat  transfer  could  be  made  a monotonically  increasing  function 
of  the  wall  temperature,  stable  operation  would  he  obtained  at  any  heat  flux  that  requires  a surface 
temperature  below  the  melting  point  of  the  wall.  Thip  condition  (Cf.  Carve  C cf  Fig.  8)  could  be 
realized  if  some  of  the  eiectricai  power  supplied  to  the  wall  were  dissipated  in  a manner  described 
by  Curve  A and  the  remainder  of  the  heat  were  transferred  to  the  test  fluid  illustrated  by  Curve  B. 

If  the  following  relation  i6  valid  for  the  partial-film-!>oiling  region  of  the  test  fluid. 


fAq-1 

~Aq~j 

1 - a~ 
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LATwJ 

A 
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the  addition  of  the  two  heat  fluxes  would  result  in  a curve  such  as  C;  thus  the  wall  temperature 
would  rise  only  very  little  for  a small  increase  of  electrical  power  over  the  peak  point. 

The  simplest  method  for  dissipating  ihe  excess  heat  as  required  by  Curve  A of  Figure  8 
and  in  addition  for  supplying  the  remain'!-^  heat  to  a test  fluid  is  shown  in  Figure  9.  The  type  of 
heat  transfer  shown  in  Curve  A of  Figure  8 is  obtained  by  flowing  a high-saturation-temperature 
liquid  up  through  a tube.  Fits  high-sstureties-temperaiure  allows  the  tube  wuil  tr-mrirmtu-r  ,n  re««'h 
a high  value  before  the  liquid  boils.  Thus  Curve  A would  be  a forced-convection  line  for  this  liquid. 
Because  of  its  function,  this  liquid  is  denoted  as  the  stabilizing  f.uid . The  fluid  under  study  ( test 
fluid)  flows  through  the  annulus  formed  by  the  tube  and  outer  jacket.  The  electrical  power  supplied 
to  the  tube  is  dissipated  outwaid  to  the  test  fluid  and  inward  to  the  stabilizing  fluid.  By  controlling 
the  press-^-c  c!  the  stabilizing  tluid,  it  is  passible  to  study  the  partial-film  and  complete-film  boiling 
regions  for  the  test  fluid  over  a wide  range  of  variables. 

It  was  thus  possible  with  the  use  of  this  sonny-stun  to  study  the  three  boiling  regions  in 
detail  up  to  the  melting  point  or  yield  point  of  the  metal  wall  without  danger  of  tube  failure.  This 
equipment  was  constructed  and  studies  with  distilled  water  have  been  completed  at  several  different 
conditions  of  bulk  temperature,  pressure,  and  velocity.  These  studies  were  made  in  order  to  broaden 
the  basic  understanding  or  the  boiling  phenomena  as  well  as  to  obiain  design  uuia. 


A. 


II.  APPARATUS 

Hydraulic  Equipment 

A schematic  diagram  of  the  recirculating  system  is  shown  111  Figure  10.  An  Ingcrsoll 
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Rand,  150-foot  head,  15-gai/min  pump,  powered  by  a Louis  Allis  Co.  3-phase,  60-cyc/sec,  5-horse- 
power  molar,  was  used  to  circulate  the  water.  After  leaving  the  accumulator  which  damped  the  pump 
nMamirr  naci Rations,  the  water  passed  through  a 6-foot-long,  1 -inch-diameter  pipe  that  was  electri- 
cally heated.  This  preheater  was  electrically  insulated  from  the  main  pipe  system  wiiu  unucsiuo 
gaskets.  The  power  was  supplied  to  this  heater  by  a 100-kva,  44G-v  stepdown  transformer  whose 
secondary  voltage  could  be  varied  from  20  to  40  v in  fifty-six  steps  by  means  of  knife  switches 
tapped  from  the  primary  winding  (Cf.  Fig.  11). 

From  the  preheater  the  water  flowed  into  the  test  section  and  then  to  a he.^t  exchanger  where 
it  was  cooled  to  a temperature  suitable  for  re-entry  into  the  pump.  The  cooling  tap  water  flowed 
through  an  annulus  formed  by  the  IS-foct-long,  ^-inch-diameter  test  fluid  pipe  and  u i;4-iuch- 
diameter  jacket. 

The  flow  rate  of  the  distilled  water  test  fluid  was  determined  by  using  a double  orifice 
type  of  flowmeter  and  a 200-inch  Barton  differential  pressure  recorder. 

A degassing  tank  was  installed  to  lower  the  quantity  of  gas  in  the  water  since  it  was 
known  from  Reference  2 turn  large  amounts  of  air  would  influence  the  boiling.  The  dashed  arrows 
in  Figure  10  show  the  path  for  recirculating  the  water  during  degassing.  The  aspirator  decreased 
the  pressure  to  3 inches  Ilg  while  the  storage  tank  heater  was  maintaining  the  water  at  saturation 
temperature.  After  54  hour  of  recirculating,  the  gas  content  of  the  water  was  decreased  to  a value 
below  0.3  cc/liter  as  determined  by  the  chemistry  section  at  the  jet  Propulsion  Laboratory  using 
the  Winkler  technique. 

S.  Test  Section 

These  studies  were  made  with  distilled  degassed  water  flowing  upward  in  a vertical 
annulus  formed  by  a 54-inch-diameter,  stainless  steei  heating  tube  and  a 2/4-inch-diameter  Pyrex 
jacket  (Cf.  Fig,  12).  Several  tests  were  also  made  using  a 0.61 -inch-diameter  Pyrex  jacket  (Cf. 

Fig.  13). 

A 1 -foe '.-length  of  type  347  stainless-steel  tubing  with  » C.GOC-inch-thick  wall  separated 
the  test  fluid  from  the  stabilizing  fluid.  Two  cylindrical  copper  bus  bars  were  slipped  over  »*•. 
tubing  and  silver  soldered  3 inches  apart.  This  spacing  provided  the  3-inch  heated  length  of  the 
test  section.  The  test  section  heating  tube  and  assembly  details  are  shown  in  Figure  14.  The 
stabilizing  fluid  was  sealed  at  the  entrance  and  exit  to  the  tube  with  robber  O-rings.  This  arrange- 
ment was  adequate  to  prevent  the  high  pressure  stabilizing  fluid  from  leaking  and  permitted  the  tube 
!o  slid*  in  the  O-rings  to  Drovide  tor  thermal  expansion. 

The  low  value  of  the  ratio  of  tube  heated  length  to  diameter  was  necessary  because  ot  the 
electrical  power  limitations  and  also  the  ciwire  to  prevent  large  changes  in  liquid  bulk  temperature 
between  the  inlet  and  outlet  of  the  test  section.  Calibration  teste  irdic“t»d  ,n  annreriable  change 
in  the  heat  transfer  coefficient  along  the  heating  tube  (Cf.  Section  IVA.). 

C.  Electrical  Power  Supply 

it  was  found  that  a 60-cyc/sec  ac  power  supply  would  cause  a 120-cyc/sec  growth  and 
collapse  of  the  bubbles,  Uonsequestiy  the  power  suppl/  chosen  for  ibis  study  was  obtained  from 
two  dc  welding  generators  (manufactured  by  the  Hobart  Co.)  which  generated  dc  power  with  the 
amplitude  of  the  voltage  ripple  less  than  one  per  cent  from  no  load  to  rated  capacity.  Regulation 
of  the  field  current  with  a rheostat  provided  smooth  control  from  no  load  to  rated  capacity  of  24 
kw  for  each  generator. 

D.  Thermocouples 

The  temperature  rise  of  both  the  tost  fluid  and  the  stabilizing  fluid  was  determined  from 
two  six-junction  copper  constantan  difference  thermocouples.  These  thermocouples  were  insulated 
and  then  inserted  into  %-incu  stainless  steel  tubes.  The  tubes  were  placed  along  the  cents?  line 
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of  tiic  entrance  and  exit  pipes  with  u 3-inch  length  irrmersed  in  an  isothermal  region.  A mixing  chamber 
was  used  downstream  from  ihe  test  section  so  thnt  the  thermocouples  read  a true  mixed  meuit  temper- 
ature (CG  Fig.  15). 

In  addition  to  the  two  difference  thermocouples,  two  copper  con6tantan  thermocouples  were 
used  to  measure  the  islet  temperature  of  the  stahili  zing  and  test  ^lutos.  .^*11  tttermo c o up les  ivsi s c*«» * 
hraled  at  the  ice  point  and  boiling  point  of  water.  While  immersed  in  wet! -stirred  baths  the  difference 
thermocouples  were  checked  over  a range  of  temperatures  against  standard  thermometers. 

An  ice  bath  was  used  as  a cold-junction  reference,  and  all  thermocouples  were  read  on  a 
manually  balanced  Rubicon  potentiometer. 


E.  Photographic  Technique 

Motion  pictures  of  the  tliree  types  of  boiling  were  taken  through  the.  Pyrex  junket  at  3wu 
frames,  sec  using  n 16-mm  Fastax  enmera.  The  bubbles  were  illuminnted  by  a General  Electric  B-I16 
Mercury  vapor  lamp  (Cf.  Fig.  16).  A concnve  mirror  placed  behind  the  lamp  focused  the  light  falling 
on  the  tube  into  a strip  % inch  wide  and  2 inches  long.  This  arrangement  provided  sufficient  light 
intensity  to  give  proper  exposure  on  F.a«tn«n  Kodak  Super  XX  safety  film  with  the  camera  lens  stop 
set  below  f/11. 


III. 


TEST  PROCEDURE 


A.  Wall-Temperature  Measurement 

The  wall  temperature  had  to  be  determined  without  disturbinK  either  the  test  fluid  or  the 
stabilizing  fluid.  A calibration  was  obtained  by  flowing  the  stabilizing  fluid  through  the  vertical 
tube  without  ihe  surrounding  test  fluid.  Three  copper  constantan  thermocouples,  electrically 
insulated  from  the  dc  voltage  gradient  in  the  tube  by  using  a 0.0005-inch-thick  mica  sheet,  were 
wrapped  around  the  outside  of  the  tube  and  spaced  one  inch  apart.  Natural  Sil-O-Cel  insulation 
was  placed  over  the  thermocouples,  and  the  entire  assembly  wns  clamped  to  the  tube  (Cf.  Fig.  17). 
This  arrangement  measured  wall  temperatures  within  4 per  cent.  To  check  the  conduction  errors  of 
the  readings,  hot  water  was  passed  through  the  tube  with  the  electrical  power  off.  The  wall  thermo- 
rc  :ples  agreed  within  2 per  cent  with  thermocouples  immersed  in  the  water  which  was  as  hot  as 
4'Jir i*  • when  u cnccK  oi  me  cicctricoi  iosuiuiin^  property  ox  me  construction  v*ss  rn due  by  fcvcrSiug 
the  power  leads,  no  shift  was  noted  in  the  thermocouple  readings. 

A series  of  tests  wns  nerformed  to  determine  the  forced  convection  line  for  the  oiokiliTm. 

O 

fluid.  The  data  obtained  agreed  within  10  per  cent  with  the  predirted  values  using  the  Sieder-Tate 
equation. 

.After  the  forced  convection  curve  was  obtained  for  the  stabilizing  fluid,  the  three  wall 
thermocouples  were  removed.  !t  was  t.hen  possible  to  determine  the  wall  temperature  without  these 
thermocouples  by  entering  the  forced-convection  curve  with  the  known  heat  flux  «s  determined  bv 
the  temperature  rise  indicated  by  the  stabilizing  fluid  difference  thermocouples.  To  obtain  the 
temperature  of  the  wall  exposed  to  the  test  fluid,  it  was  necessary  to  correct  for  the  temperature 
change  through  the  wall.  The  position  of  zero  temperature  gradient  was  obtained  by  assuming 
constant  electrical  resistance  through  the  wall.  Then  by  Figure  18 
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Assuming  constant  electrical  and  thermal  resistance,  the  temperature  drop  from  J.\  to  the  teat  fluid 
interface  is  given  Ly  McAdams  (Cf.  Hcf.  2)' 
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The  temperature  variation  in  the  wall  as  calculated  by  relation  (3)  or  (4)  was  always  less  than  45° F; 
thus  the  assumption  of  constant  electrical  and  thermal  resistance  was  reasonable. 


B.  Stabilizing  Fluids 


Distilled  water  at  a pressure  of  4b0  psia  and  flowing  at  40  ft/acc  was  used  as  the  stabilizing 
fluid  for  wall  temperatin-es  up  to  500“ F.  This  upper  limit  resulted  when  the  water  began  to  boil  isad 
thus  maintained  the  wall  temperature  substantially  constant  regardless  of  any  additional  electrical 
power  input.  Sirce  the  wall  temperature  remained  a fixed  amount  above  the  boiling  pcin;  of  the 
stabilizing  fluid,  higher  wall  temperatures  could  have  been  reached  by  increasing  the  pressure  of 
the  stabilizing  -“nter.  However,  the  small  change  in  saturation  temperature  with  pressure  made  this 
tol  jlion  impractical.  As  an  alternative,  nroiher  higher-- boiling-point  liquid  or  a high-pressure  gas 
could  be  used  instead  of  the  water  as  the  stabilizing  fluid.  Dry  nitrogen  gas  proved  to  be  a satis- 
factory stabilizing  fluid. 

Ten  1.5-cu  ft  tanks,  initially  filled  with  2200-psia  dry  nitrogen  gas,  were  connected  to  the 
system.  A pressure  regulator  delivered  500-psia  nitrogen  to  the  heater  tube.  After  passing  through 
the  test  section,  the  nitrogen  was  discharged  at  sonic  velocity  to  the  atmosphere. 


A potential  difficulty  with  the  use  of  nitrogen  gas  as  the  stabilizer  fluid  was  the  isenthalpic 
ehange  in  temperature  with  pressure.  It  was  possible,  however,  to  employ  the  nitrogen  gau  since  the 
stagnation  temperature  of  the  gas  entering  the  t.cmt  section  verte-’  ih««  TV'C  Ju. Ir.a:  the  test. 
Reference  to  F'gure  19  indicates  that,  b>  ;~g  -wing  the  u.^-ogen  pressure  at  500  psia,  the  discharge 
temperature  hi*  eimerted  to  remain  between  -t  2°C  and  — 32°C  as  the  storage  tank  pressure 

decreased  from  2200  psia  to  1000  psia. 


The  forced  convection  cirvt  for  the  nitrogen  gas  was  determined  in  the  same  manner  as  for 
the  water  stabilizing  fluid.  For  a constant  power  input  to  the  rube,  the  see- sure  d wall  temperature 
decreased  less  than  10  per  cent  ae  the  storage  tank  pressure  decreased  from  2000  psia  to  1000  psia. 
As  a consequence,  the  heat  transfer  was  assumed  in*.iepv  t-ie-t  of  storage  tank  pressure  when  the 
pressure  in  the  tank  was  above  1000  psia. 

The  forced  convection  curve  was  slightly  concave  downward  because  of  the  higher  exit 
temperature  of  the  nitrogen  from  the  test  section  at  higher  heat  flux.  The  higher  exit  temperature 
resulted  in  <*  lower  density  [p~(l/T)l.  The  exit  velocity  remained  sonic  hut  the  velocity  of  sound 
increased  with  temperature  (c^^/T).  The  net  effect  ivcs  a lower  mass  flow  rate  [G~  (l/-,/T).  The 
lower  mass  flow  caused  a decrease  in  the  gaa  film  coefficient  and  resulted  in  a concave  downward 
forced  cunvcctiGu  curve* 


A power  balance  was  obtained  by  measuring  the  temperature  rise  through  the  test  section 
of  the  test  fluid  and  the  water  stabilizing  fluid.  The  heat  transfer  calculated  from  these  two  values 
agreed  with  power  measurements  within  five  per  cent.  Because  of  the  difficulty  in  measuring  th-"> 
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temperature  rise  of  the  nitrogen,  no  attempt  was  made  to  obtain  a power  balance  when  the  nitrogen 
gas  was  used.  The  heat  transferred  to  the  nitrogen  was  obtained  by  subtracting  from  the  total  power 
input  the  heat  absorbed  by  the  test  water.  Th~  power  balance  for  ‘he  previous  tests  with  water  as 
the  test  and  stabilizing  fluid  indicated  that  this  procedure  should  give  heat  flux  values  accurate  to 
five  per  cent. 

£.  Power  Measurements 

Thi  power  delivered  to  the  test  section  was  measured  with  an  ammeter  and  voltueter,  each 
accu  ite  to  ne  per  cent  of  full  scale  reading.  For  convenience,  the  voltage  measurement  wrb  aict-n 
6-1  'he  outer  ends  of  the  copper  bus  bars  instead  of  across  the  heated  length  of  the  tub  This 
arrangement  resulted  in  less  than  a one-per  cent  error-  The  total  estimated  error  in  power  measurement 
wae  less  than  3 per  cent.  The  heat  flux  was  calculated  from  the  measured  power  value  and  tube-area 
measurements.  The  heat  flux  was  known  within  an  uncertainty  of  <4  per  cent.  This  method  was  checked 
against  values  obtained  by  measuring  current  flow  and  calculating  power  input  from  the  product  I2R. 
Because  of  a 3-per  cent  variation  in  tube  wall  thickness  between  various  tubes,  the  second  method 
was  discarded  in  favor  of  measuring  both  current  and  voltage. 

IV,  DISCUSSION  OF  RESULTS 

A.  Th*  Mechanism  of  Transition  Between  the  Three  Boiling  Regions 

In  an  effort  to  understand  the  mechanism  by  which  heat  is  transferred  from  a hot  metal  wall 
to  a subcooled  boiling  liquid,*  high  speed  motion  pictures  were  taken  of  the  three  boiling  regions. 
Emphasis  was  placed  or.  the  partial  and  complete  film  boiling  mechanisms;  die  second  part  of  this 
report  on  boiling  beat  transfer  will  deal  with  the  complex  nucleate  boiling  process. 

A brief  discussion  of  the  hailing  mechanism  is  gives  is  the  introduction  of  this  paper,  but 
the  details  of  the  mechanism  by  which  the  boiling  process  changes  from  nucleate  to  partial  film 
and  finally  to  complete  fiim  boiling  are  now  discussed. 

The  first  problem  to  be  encountered  in  understanding  ttie  boiling  mechanism  is  the  process 
by  which  a bubble  is  generated.  Since  this  problem  is  beyond  the  scope  «»*  she  present  study,  the 
commonly  accepted  conditions  are  assumed  to  the  effect  that  the  bubbles  may  form  ouiy  at  nuclei 
of  sufficient  size  to  overcome  the  initial  surface  tension  force  of  the  liquid  end  that  these  nuclei 
are  made  up  of  air  or  vapor  or  a combination  of  both.  The  size  of  nuclei  which  is  required  to  initiate 
boiling  in  any  liquid  under  given  conditions  decreases  with  increasing  liquid  temperature,  T'hu«  for 
any  size  distribution  of  nuclei  in  a liquid,  there  is  a greater  number  of  active  nuclei  at  higher  liquid 
temperatures  and  consequently  a greater  number  of  bubbles.  The  following  discussion  is  based  on 
the  premise  that  the  population  of  bubbles  inc.'cases  with  the  temperature  of  the  liquid.  A more 
detailed  discussion  of  the  generation  of  bubbles  is  given  in  Port  II,  Section  IV. 

For  simplicity  in  discussing  the  boiling  premess,  it  is  considered  that  a vertical  metal 
wall  is  immersed  in  a stagnant  poo!  of  degassed  subcooled  water  at  1 -atmosphere  pressure.  Since 
not  all  of  the  dissolved  gas  has  bee:;  removed  from  the  water,  small  bubbles  (fow  near  the  hot 
wall  from  the  nuclei  that  ere  sufficiently  large  to  overcome  the  lirruid  surface  tension.  The  first 
noticeable  effect  of  bubbles  on  heat  transfer  appeers  in  average  degassed  water  (0.3  cc  of  air  per 
liter  of  water)  when  the  wall  temperature  exceeds  the  liquid  saturation  temperature  by  about  30°F. 
These  ujcfaate  boiling  bubbles  grow  and  collapse,  thereby  creating  turbulence  in  the  liquid  near 
the  wall.  i he  hot  liquid  near  the  wail  is  iorced  away  by  this  turbulence  and  replaced  by  the  cooler 
bulk  liquid.  The  temperature  of  the  local  eddies  of  water  may  then  alternate  from  superheated  to 
below  saturation-  Since  the  wall  ienmerature  is  only  slightly  above  the  saturation  temperature,  a 

•The  te>.n  subcoolsd  boiling  liquid  is  con».or,!y  used  for  simplicity  despite  the  obvious  ambiguity. 

The  term  replaces  ihe  longer  statement  “the  liquid  whose  built  i«:im>«-.r»t>ire  is  below  saturation  temperature 
'out  whose  boundary  layer  is  superheated  so  that  bubbles  may  iorre  next  to  the  heat  transfer  surface". 
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a small  uumber  of  widely  spaced  bubbles  may  maintain  the  liquid  near  tbe  wall  between  the  bubbles 
»t  h temperature  close  to  saturation.  Whenever  the  liquid  temperature  increases  to  the  value  needed 
to  evaporate  iuto  a nucleus  which  is  present,  a new  btihbie  forms. 

As  the  wall  temperature  is  increased,  the  bubble  population  al  Sv  increases  since  more 
nuclei  ere  of  the  size  required  to  initiate  a bubble.  The  population  of  bubbles  increases  to  the  number 
that  is  needed  to  create  sufficient  turbulence  to  cool  the  remaining  liquid  nexi  to  the  wall  be  lev  it" 
nucleation  temperature.  When  the  wall  is  about  10(FF  above  the  water  saturation  temperature,  the 
population  has  increased  so  that  at  maximum  size  the  bubbles  cover  100  per  cent  of  localized  areas 
of  the  hot  wall.  These  localized  areas  are  entirely  covered  by  about  a dozen  bubbles  which  coalesce 
to  form  a locul  vapor  blanket.  At  any  in3tant  these  local  areas  that  are  entirely  covered  with  bubbles 
may  occupy  only  30  per  cent  of  the  total  area.  At  the  peak  heat  flux  point,  about  30  per  cent  of  the 
totn!  wall  area  is  occupied  by  vapor  in  the  form  of  blankets  and  individual  bubbles  (Cf.  Fig.  21). 

Figure  22  illustrates  the  vertical  wall  with  a local  film  blanket  for  tbe  condition  al  the 
transition  from  nucleate  tc  partial  film  boiling.  The  buoyant  force  on  this  large  vapor  blanket  exceeds 
the  adhesion  force  and  tend),  tc  slide  the  blanket  upward  along  the  vertical  wall.  However,  since  the 
stagnant  liquid  near  the  wall  above  the  blanket  resists  motion  more  strongly  than  tbe  liquid  far  from 
the  wall,  the  vapor  will  float  out  into  the  liquid  as  shown.  When  the  vapor  detaches  from  the  wall, 
cool  liquid  is  drawn  in  as  replacement. 

These  thick  local  blankets  of  vapor  alternately  formfrom  coalescing  nucleate  bubbles  and 
float  out  into  the  main  body  of  the  liquid  to  condense.  The  partial  film  boiling  continues  until  the 
wall  temperature  reaches  about  1000° F above  the  saturation  temperature..  At  this  high  temperature, 
the  bubble  population  has  increased  to  the  point  where  the  bubbles  contact  one  another  when  only 
partially  grown.  For  this  condition  the  many  local  film  blankets  now  occupy  ICO  per  rent  of  the 
entire  wall  and  coalesce  tc  form  a complete  film  (Cf.  Fir.  21).  Since  the  entire  wall  is  covered,  the 
vapor  is  forced  upward  by  the  buoyant  forces  rather  thar.  floating  away  from  the  wall  (Cf.  Fig.  23). 

The  film  thickness  is  now  adjusted  so  that  the  vapor  flowing  upward  is  just  replaced  by  the 
evaporation  into  the  film. 

It  is  seen  that  the  individual  nucleate-boiling  type  of  bubble  plays  a prominent  role  in  nil 
three  boiling  regions.  The  only  difference  between  the  various  boiling  regions  is  the  population  of 
nucleate  bubbles  a few  milliseconds  after  the  wall  is  wetted. 

A discussion  of  the  reverse  process  of  starting  with  a hot  wall,  which  is  cooled  down  to 
the  ronbotitng  temperature.  .uu,  «id  la  usd'rstsnding  the  mechanism.  This  process  would  be  the 
case  of  quenching  a hot  metal  body  in  a subcooled  liquid. 

At  the  instant  when  the  hot  body  (at  least  1000'F  ubove  the  liquid  saturation  temperature) 
is  immersed  in  the  subcooled  liquid,  nucleate-boiling  bubbles  form.  After  a cycle  or  wo  of  growth 
and  collapse,  the  bubble  population  is  so  greet.  (Ci.  Fig.  24)  that  several  bubb.es  coalesce  into 
local  blankets  of  vapor  which  in  turn  coalesce  into  a complete  vapor  film.  The  liquid  evaporates 
into  the  film,  and  the  vapor  flows  upward  by  buoyant  forces.  The  film  thickness  adjusts  so  that 
the  evaporation  into  the  film  just  equals  the  mass  of  vapor  flow  out  of  the  film.  Heat  is  conducted 
and  radiated  from  the  wall  through  the  vapor  film  to  the  liquid.  As  the  wall  is  cooled,  the  vapor- 
film  thickness  must  decrease  to  maintain  the  balance  between  evaporation  into  the  film  and  mass 
flow  out  of  the  film.  The  heat  transfer  by  conduction  through  the  film  decreases  since  the  temper- 
ature difference  Tw-TBal  decreases  faster  than  the  film  thickness  (Cf.  Eq,  A-10).  !r.  addition,  the 
heat  transfer  by  radiation  decreases.  The  totcl  heat  transfer  must  consequently  decrease  as  the 
temperature  of  the  wall  is  lowered. 

The  film  bailing  in  the  pool  has  caused  turbulence  in  the  water.  These  water  currents  may 
at  times  penetrate  the  vapor  film  and  wet  the  wall.  If  the  waL'  is  sufficiently  hot,  nucleui::  bubbles 
form  over  100  per  ceui  cf  the  wetted  area,  and  the  complete  film  is  again  established  (Cf.  Fig.  24). 
However,  if  the  film  is  penetrated  by  liquid  when  the  wail  t.a  no  longer  hot  enough  to  produce  100 
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per  cent  coverage  of  nucleate  bubbles  over  the  entire  wetted  urea,  the  resultiug  partial  film  boiling 
creates  additional  turbulence  which  breaks  up  more  of  the  film.  The  film  may  be  broken  at  any  place 
during  the  transition  from  complete  io  partial  film  Lolling  although  the  liquid  is  ny'rc  likely  to  wet 
at  the  bottom  of  the  wall  where  the  film  is  thinnest.  The  liquid  at  this  point  then  boils  and  creates 
turbulence  which  forces  the  cool  liquid  against  the  remaining  film.  This  action  causes  more  of  the 
well  in  he  wetted  and  more  nucleate  bubbles  to  form.  These  nucleate  bubbles  have  a population 
that  allows  local  blankets  of  vapor  to  be  formed  thus  cr»»*ing  more  turbulence,  in  this  manner  the 
complete  film  is  removed  from  the  wail  and  replaced  with  local  blankets  in  the  partiai*fiim*boi!ing 
region.  Heat  continues  to  flow  from  the  wall  into  the  subcooled  liquid  by  the  mechanism  of  induced 
turbulence  and  latent  heat  transport  us  the  blankets  float  away  from  the  wall. 

The  vtal!  continues  t«  be  cooled,  and  the  population  of  lccal  vapor  blankets  decreases, 
thus  allowing  more  time  for  the  intermediate  nucleate  bubbles.  This  decrease  in  the  population  of 
vapor  blankets  results  in  a higher  heat  transier.  When  the  wall  is  couied  t?  the  point  where  the 
bubble  population  is  too  small  to  allow  coalescence,  beat  is  transferred  entirely  by  the  nucleate- 
boiling  mechanism.  As  the  wall  temperatire  is  decreased  further,  the  nucleate-bubble  population 
decreases,  and  the  heat  transfer  consequently  decreases.  Finally  the  wall  temperature  fails  below 
the  liquid  saturation  temperatiye.  znd  the  body  is  cooled  to  the  liquid  temperature  by  fret 
convection. 

In  addition  to  these  conclusions  regarding  the  mechanism  of  boiling,  the  following  experi- 
mental results  have  been  obtained. 

B.  Forced  Convection 

The  heat-transfer  curves  obtained  for  the  forced-convection  nonboiling  region  (Cf.  Figs. 

25  through  28)  agreed  within  10  per  cent  with  values  calculated  from  the  Siedcr-Tate  equation. 


i,C^Re0-8  ( TW ~ TL ) 

«*  t/  \ r-iMi  / 


where  the  best  value  of  C was  0.025  (Ct.  Fig.  29). 

C.  Huclsato  Boiling 

The  nucleate-boiling  data  correlated  in  Figure  30  were  vyithin  5 per  cent  of  the  values 
ffivFr.  nv  McAdams  (Cf.  Ref.  2). 

The  heat  flux  in  the  nucleate-boiling  region  varied  from  incipient  boiling  to  the  peak 
heat  flux  with  substantially  a constant  wall  temperature.  Because  of  the  small  variation  of  wall 
temperature  with  heat  flux  in  this  region,  oely  two  design  farts  must  be  known  about  the  nucleate- 
boiling  process:  (1)  tbe  maximum  neot  flux  that  may  be  transferred  at  a reasonably  low  wall 
temperature  (usually  the  heat  flux  at  the  transition  from  nucleate  to  partial  film  boiling)  and  (2) 
the  wall  temperature  required  to  cause  incipient  boiling  in  the  Luuuu&ry  layer  of  the  coolant 
liquid. 

The  present  study  has  furnished  information  on  the  mechanism  of  nucleate  boiling, 
partial  film  boiling,  and  complete  film  boiling.  In  order  to  obtain  additional  information  on  the 
factors  controlling  the  maximum  heat  flux,  more  detailed  pool-boiling  studies  of  nucleate  boiling 
were  made  fur  water,  water-Aerosoi  solution,  and  carbon  tetrachloride  containing  various  amounts 
of  gas.  The  combined  program  has  resulted  in  a better  understanding  of  the  boiling  mechanism 
which  has  allowed  the  formulation  of  a semi-empirical  theory  that  successfully  predicts  the  peak 
heat  transfer  rates  (Ci.  Part  I!  Sect.  V).  The  calculation  of  peak  heat  flux  by  this  method  still 
requires  an  experimental  deteirr.isuticn  of  bubble  growth  rates. 

The  second  important  fact,  that  of  knowing  the  wall  temperature  at  which  the  liquid 
will  nucleate  boil,  is  equally  difficult  to  study.  To  data,  so  theory  is  avails ble  whit  a cciu  predict 
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successfully  the  wall  temperature  required  to  initiate  boiling.  Several  unknown  facts,  such  as 
the  average  size  and  distribution  of  gas  nuclei  or  other  foreign  suclcaiicn  part’d***;  prohibit 
an  analytical  solution  of  the  problem  at  this  time.  However,  a rough  estimate  of  wall  temperature 
may  be  made  since  experirae ntal  results  have  indicated  that  all  moderately  degassed  liquids 
tested  boil  when  the  wall  temperature  exceeds  the  liquid  bubl ! e-point  temperature  by  less  than 

100°F. 

D.  Partial  Film  Bolling 


Reference  to  the  Figures  25,  26,  and  27  shows  that  the  wall  temperature  in  the  partial- 
film-boiling  region  is  effected  by  pressure  and  subcooling  in  al  out  the  samf  manner  as  in 
nucleate  boiling.  The  wall  temperature  increases  with  pressure  so  that  the  difference  between 
wall  temperature  and  saturation  temperature  remains  approximately  constant  at  any  given  heat 
flux.  This  correlation  is  shown  in  Figure  31,  where  heat  flux  is  plotted  against  the  difference 
between  wall  temperature  and  liquid  saturation  temperature. 

The  high-speed  motion  pictures  of  this  type  of  boiling  indicate  that  the  cycle  consists 
of  the  following  events:  (1)  nucleate  bubbles  form;  (2)  the  population  of  the  bubbles  increases 
until  several  coalesce  into  local  blank****  of  vapor;  (3)  these  blankets  grow  out  into  the  cooler 
liquid  away  from  the  hot  wail  auu  are  then  condensed  while  still  attached  to  the  wall  or  float 
away  by  buoyant  fores  and  condense  in  the  main  stream  (Cf.  Fig.  22). 


E.  Fiim  Boiling 

Film  boiling  may  occur  in  several  different  types  of  equipment  such  as:  nuclear  reactors, 
condensers,  binary  me  rcury-w  Liter  boilers,  and  quenching  apparatus. 

The  heat  transfer  film  boiling  was  found  to  he  independent  of  the  water  pressure, 
velocity,  or  subcooling  over  rue  range  investigated  (Cf.  Figs.  25  through  28). 

The  water  flow  at  5 ft/sec  was  reversed  in  one  test  so  as  to  flow  from  top  to  bottom. 

The  vapor  film  continued  to  flow  upward.  This  result  indicated  that,  ior  a vertical  heat  transfer 
surface  and  the  range  of  velocities  investigated,  the  main  driving  force  for  the  vapor  film  was 
huoyarice,  not  pressure  drop  due  to  friction  cr  shear  stress  of  the  water  flow.  A correlation  for 
the  complete-film-boiliog  heat  transfer  from  a vertical  wall  is  presented  in  Figure  32a. 

Measurements  of  the  photographs  indicated  that  the  vapor  velocity  was  about  25  ft/oec 
regmuiess  oi  the  varir-J  €•••?*' t linns.  This  v-lsse  H-rrees  with  the  calculated  estimates  given  in 
the  Appendix. 

!n  th*  investigation  of  film  boiling,  the  electrical  power  supplied  to  the  heating  tubs 
was  increased  until  the  tube  failed.  F'ailure  was  not  caused  by  melting  of  the  entile  tube-  hut 
usually  resulted  from  the  decrease  in  allowable  stress  at  high  wall  temperature.  T be  lower 
yield  point  could  not  support  the  internal  pressure,  and  the  tube  bulged  as  shows!  ia  Figure  33n. 
Another  less  frequent  cause  of  failtre  was  a local  burnout  of  about  yu-inch  diameter  caused  by 
an  imperfection  in  the  tube  wall.  In  such  a case  the  stabilizing  fluid  flowed  out  of  the  pinhole 
and  prevented  the  burnout  from  spreading  (Cf.  Fig.  33a). 


F.  Accuracy  of  Data 

The  data  which  were  obtained  from  this  study  were  probably  within  10  per  cent  of  the 
true  hv?at  flux  values  in  the  forced-convoctioa,  nucleate-boiling,  and  complete-film-Loiling 
region.  The  accuracy  decreased  to  about  15  per  cent  in  the  p«rtinl-film-bei!iag  region  because 
of  the  difficulty  of  maintaining  a constant-mass-flow  rate  with  the  violent  growtn  and  collapse 
of  the  large  vapor  rja  uses. 

V.  CONCLUSIONS 


The  major  aiin  of  this  part  of  the  investigation  was  to  study  the  three  boiling  regions 
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in  order  to  understand  better  the  mechanism  of  boiling  heat  transfer.  This  program  has  yielded 
sutticient  information  to  indicate  how  the  boiiing  process  changes  from  nucleate  into  partial  film 
boiling  and  finally  iuto  complete  film  boiling  as  the  temperature  of  the  heat-transfer  surface  is 
raised.  The  peak  heat  flux  values  were  seen  to  occur  approximately  =;  the  transition  from  nucleate 

iO  partial  liiui  DOillu^t 

The  equipment  that  was  developed  for  this  program  proved  satisfactory  for  studying  the 
"unstable**  partial-film-boiling  region  as  well  as  the  complete-film-fcoiling  region.  With  slight 
modification  to  allow  for  the  use  of  a higher-pressure  stabilizing  gas,  the  apparatus  may  be  used 
to  study  the  boiling  regions  of  water  flowing  at  hippier  velocities  than  were  employed  in  ihia 
program  (e.g.,  the  use  of  1000-psia  nitrogen  gee  would  stabilize  water  flowing  of  any  velocity  up 
to  25  ft/sec).  Other  liquids  that  have  lower  peak  heat-flux  values  than  water  may  be  stabilised 
with  nitrogcu  gas  at  much  lower  pressures. 

PART  II 

A DETAILED  PHOTOGRAPHIC  STUDY  OF  NUCLEATE  POOL  BOILING 

SUMMARY 

The  second  phase  af  the  investigation  consisted  of  s study  of  poo!  boiling  using  distilled 
water  and  commercially  pure  carbon  tetrachloride.  The  experimental  program  was  aimed  at  ohtaining 
fundamental  information  on  the  behavior  of  nucleate  bubbles  forming  on  a stainless-steel  heating 
strip  and  the  role  they  play  in  boiling  heat  transfer.  High-speed  motion  pictures  were  taken  of  the 
nucleate  bubbles  at  a liquid  pressure  of  one  atmosphere,  liquid  temperature  range  from  170°F 
helow  saturation  to  saturation  temperature,  and  heat  fluxes  from  incipient  boiling  to  the  transition 
from  nucleate  to  partial  film  boiling.  In  addition,  the  effects  of  dissolved  gas  »nd  «t»rf«re  tension 
on  kuhhle  dynamics  were  studied. 

As  a result  of  this  program  it  has  been  possible  to  propose  r meciianiBiE  for  the  growth 
and  collapse  of  nucleate  buhhles.  An  empirical  expression  has  heen  obtained  that  relates  the 
measured  huhble  velocity  and  size  to  the  heat  flux  at  the  transition  from  nucleate  to  partial  film 
boiling.  This  relation  is 

Nu  = 0.053  Re°-8Pr 

The  velocity  and  effective  diameter  appearing  in  this  equation  still  have  to  he  determined  ex- 
perimentally; all  other  factors  are  liquid  properties.  The  values  of  the  peak  heat  u-ouafer  caiuu- 
lated  from  this  relation  agree  within  15  per  cent  with  the  experimental  data. 
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The  high  combustion  gas  temperatures  nnd  gae  film  beat  transfer  coefficients  in  rocket 
motors  result  in  heat  transfer  rates  far  in  excess  of  the  values  encountered  in  the  conventional 
industrial  heat  transfer  equipment.  Heat  fluxes  of  the  order  of  10  Btu/sq  in. -sec  are  not  uncommon 
in  the  modern  rocket  motors.  To  cool  these  motors  hy  regenerative  forced  convection  would  require 
excessively  high  velocities  with  n resulting  high  pressure  drop.  In  addition  the  coolant  passages 
would  in  many  cases  be  so  small  that  maintaining  the  required  tolerances  would  cause  undesirable 
m«n>jfnc?'_r.:,!g  difficulties.  It  is  possible,  however,  to  transfer  these  high  heat  fluxes  hy  nucleate 
boiling  at  much  lower  velocities.  Reference  to  Figure  26  shows  that  for  the  given  conditions  ten 
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bailing  forced  convection  at  the  same  free  stream  velocity.  The  wall  temperature  increased  only 
i(f  F to  result  in  this  tenfold  increase  in  hert  flux.  Since  this  small  temperature  increase  could 
not  change  the  physical  properties  of  the  rocket-motor  wall  to  any  extent,  the  only  important 
design  dots  needed  are  the  wall  temperature  required  tc  initiate  nucleate  boiliug  and  the  peak  heat 
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flux  which  may  be  transferred-  Although  the  variation  of  heat  flux  with  wall  temperature  in  the 
nucleate  boiling  region  is  of  little  importance  for  design  purpose?,  *t  >®  of  interest  in  understanding 
the  process  of  boiling  heat  transfer  as  will  be  seen  in  the  discussion  of  bubble  dynamics. 

In  addition  to  the  application  of  nucleate  boiling  to  the  cooling  of  rocket  motors,  the 
process  is  of  importance  in  the  design  uf  nuclear  reactors  and  the  proper  quenching  of  metals. 

Although  the  conventional  heat  exchangers  and  condensers  usually  are  not  limited  in  sire  and 
weight  to  the  same  extent  as  flight  equipment,  it  would  be  advisable  economically  to  design  this 
equipment  to  operate  in  the  nucleate-boiling  region  also  whenever  possible. 

B.  Lit 

The  more  significant  advances  in  the  field  of  pool-boiling  heat  transfer  up  to  the  present 
will  be  discussed  briefly  in  this  section  of  the  introduction. 

Mosciki  and  Broder  (Cf.  Ref.  1)  were  the  first  to  investigate  the  variation  of  heat  flux 
with  wall  temperature  in  the  free  convection  and  nucleate  boiling  region  by  employing  an  electrically 
heated  wire  in  a pool  of  water.  Nikiyama  (Cf.  Ref.  5)  extended  the  study  tor  bulk-boiling  water  to 
include  the  complete  film  boiling  region. 

Drew  and  Mueller  (Cf.  Ref.  8)  postulated  that  in  the  region  of  rapid  increases  in  heat  flux 
for  small  increases  in  wall  temperature  the  liquid  boiled  on  the  heating  surface  at  aciive  nuclei. 

This  proposal  resulted  in  terming  the  boiling  region  nucleate  boiling. 

McAdams,  et  al,  (Cf.  Rtf.  9)  found  that  the  beat  transfer  in  nucleate  boiling  was  independent 
of  wire  diameter  from  0.016  to  0.024  inch.  For  wire  diameters  of  0.004  to  0.008,  McAdams  stated  that 
at  any  given  beat  flux  the  temperature  difference  between  the  wall  and  the  water  increased  for  larger 
wire  size.  (It  should  be  noted  that  these  conclusions  do  not  include  the  peak  heat  flux  values.) 

Jakob  (Cf.  Ref.  7)  proposed  that  the  increase  in  heat  transfer  in  the  nucleate  boiling  region 
for  bulk  bailing  vms  due  »o  the  agitation  of  tbe  liquid  near  the  wall  caused  by  the  detaching  bubbles. 
This  conclusion  was  reached  after  showing  that  tbe  latent  heat  transport  mechanism  alone  was  in- 
sufficient to  account  for  tbe  hig^i  beat  flux  obtained  by  nucleate  boiling. 

King  (Cf.  Rsf.  10)  found  that  the  wall  temperature  exceeded  the  liquid  saturation  temperature 
by  an  amount  which  decreased  with  the  addition  of  dissolved  gases  to  tbe  liquid,  the  presence  of 
solid  impurities  in  the  liquid,  or  by  increasing  tbe  pressure  of  the  liquid. 

Sauer  (Cf.  Rsf.  11),  Akin  end  ?«c A •!«»-•>•  (Cf.  Ref.  12).  and  Farbcr  and  Scarab  (Cf.  Ref.  13) 
showed  that  different  metal  heating  surfaces  produced  appreciable  variations  in  tbe  nucleate-boiling 
heat  transfer  for  any  giyen  temperature  difference  between  the  wall  uud  liquid.  According  to  Kaniakis 
and  Sherman  (Cf.  Ref.  14)  the  surface  roughness  of  the  heating  surface  did  not  influence  the  vulu^  of 
the  peak  heat  flux. 

Gunt'ncr  (Cf.  Ref.  IS)  employed  r.  high  speed  camera  study  nucleate  boiling  and  speculated 
that  the  high  rates  of  heat  transfer  were  due  to  induced  liquid  turbulence. 

In  addition  to  t'ue  interest  in  bubble  dynamics  as  affecting  boiling  beat  transfer,  considerable 
effort  has  been  expended  in  tbe  study  of  cavitating  flow.  Unlike  tbe  preceding  studies  of  boiling  heat 
transfer,  the  hydraulic  ca'd'.ating  flow  has  been  studied  in  detail,  both  analytically  and  experimentally, 
by  many  investigators,  including  Knapp  ano  Hollander  (Cf.  Ref.  16)j  Harvey,  et  al  (Cf.  Ref.  17); 
Plesset  (Cf.  Ref.  18);  Ellis  (Cf.  Ref.  19);  Dergarabedian  (Cf.  Ref.  20);  and  Gilmore  (Cf.  Ref.  21). 

C.  Purpose  of  Present  Study 

The  study  of  nucleate  boiling  heat  t~n“fer  may  be  divided  into  three  phases: 

1.  Bubble  generation 

2.  Bubble  growth  and  COtiopS  C 

3 ? Influence  of  bubble  dynamics  on  beal  transfer 
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The  subject  may  be  subdivided  further  into  the  categories  of  beiiiag  in  a poo!  and  boiling 
with  forced  reave  crion. 


The  phases  cbosen  for  this  experf'oo*'  1 'nvestK'uiion  were  bubble  growth  and  collapse  as 
well  as  the  influence  of  bubble  dynamics  on  bent  transfer  in  pool  boiling.  Two  liquids  with  large 
differences  in  physical  properties  e chosen  for  this  investigatio:::  water  and  carbon  tetra* 
chloride,  in  addition,  since  preiiminan’  investigation  a had  indicated  iliut  the  surface  tension  and 
gas  content  of  the  liquid  played  an  important  role  in  the  nucleate  boiling  process,  a water-Aerosol 
solution  and  water  containing  different  amounts  of  air  were  investigated. 


II.  DETERMINATION  OF  PEAK  HEAT  FLUX 


Description  of  Equipment 

1.  Purpose.  The  discussion  contained  in  the  introduction  to  this  phase  of  the  report 


explained  the  desirability  of  knowing  the  maximum  heat  flux  that  could  be  transferred  to  any  liquid 
by  the  process  of  nucleate  bo’l'ng:  Tbi«  point,  which  oernre  where  the  heat  transfer  mechanism 
changes  from  nucleate  to  partial  film  boiling,  was  determined  by  increasing  the  power  supplied  to 
the  heat  transfer  strip  until  the  surfsce  melted.  The  factors  responsible  for  this  failure  of  the  heat 

8urfscs  have  been  discuss^'!  is  detail  in  Part  Is 

The  present  study  was  made  in  an  attempt  to  obtain  a basic  understanding  of  the  mechanism 
of  the  nucleate-boiling  process  that  would  aid  in  obtaining  a correlation  that  could  predict  the  peak 
heat-transfer  rates.  To  simplify  the  problem,  it  was  decided  to  confine  this  phase  of  the  study  to 
psa!  boiling.  This  condition  eliminated  the  complications  which  would  result  with  forced  convection 
where  the  nucleate  bubbles  move  downstream  at  about  8/l0  of  tbe  tree  stream  velocity  (Cf.  Ref.  4). 
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wide,  !4-inch  long  and  0.0043-inch  thick,  type  304  stainless-steel  beating  strip  immersed  in  distilled 
water.  T^e  test  f»njc  WSS  mads  of  O Pmv  dan  9fv.inrh  d!*»”>eter  rvHnder.  1 inches  lone,  with  flat 
ends  (Cf.  Fig.  32b).  The  distilled  water  was  boiled  in  a 2-liter  flask  until  one  half  of  the  water  was 
evaporated.  This  procedure  lowered  the  dissolved  gas  content  to  less  than  0.2  cc  of  air  per  liter  of 
water  as  determined  by  tbe  chemistry  section  of  the  jet  Propulsion  Laboratory  using  the  Winkler 
technique.  After  degausing,  the  water  was  drawn  by  means  of  an  aspirator  into  tbe  test  section 
where  tbe  bubbles  could  be  photographed  through  the  flat  ends  of  the  tank.  A thermometer  was  placed 
!4  inch  above  the  heating  strip  to  measure  a reference  bulk  liquid  temperature. 

The  power  supply  for  electrically  heating  the  stainless-steel  strip  was  obtained  from  a dc 
welding  generator  manufactured  by  the  Hobart  Co.  Regulation  of  the  field  current  by  the  >j»e  of  a 


rheostat  provided  smooth  control  irum  no  iuau  Lu  tuc  rt*tcu  capacity  of  24  kw  With  the  Kiupiuiiuc  of 


the  voltage  ripple  less  than  Cue  per  cent* 

The  heating  strip  was  cut  from  a roll  of  type  304  stainless  steel  with  precision  shears.  The 
"trips  chosen  for  these  tests  varied  less  than  O.OUl  inch  in  the  0.i25-incb  width  and  less  than  0.0001 
inch  :n  the  0.0043-inch  thickness.  To  determine  tbe  electrical  resistance  of  the  heating  strip,  five, 
l%°volt,  dry  cell  batteries  arranged  in  parallel  were  connected  in  series  with  a 3-ohm  resistor,  a 0.5 
per  cent  anmeter,  and  the  stainless-steel  strip.  The  voltage  drop  across  a 1-inch  length  of  the  strip 

• _ * » i .1  t 1 111  1 • . . • _ . _ _ rji|.  _ , _ . 

was  aetennmea  oy  cue  use  01  a lutnci-otmtnceu  precision  potentiometer*  i uc  avc» avc  icoiatiVity 

measured  in  this  manner  was  0*0263  ohm  per  inch  of  length*  Figure  33b  shows  tbe  curve  of  heat  flux 
v«  current  that  was  obtained  and  then  used  for  the  present  program.  The  variation  of  tbe  resistivity 
- f the  atrip  with  temperature  wea  negligibly  small  over  the  rauge  of  interest.  Tlie  heal  flux  deter- 
mined in  this  manner  was  estimated  to  be  witbin  1 per  cent  of  tbe  true  value. 

3.  Carbon  Tetrachloride.  Considerable  difficulty  was  encountered  in  an  effort  to  remove 
the  disoolved  gas  from  the  carbon  tetrachloride.  When  this  liquid  was  subjected  to  the  same  degassing 
technique  as  had  proved  satisfactory  for  the  water,  large  quantities  of  noncondensable  gas  bubbles 
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w?r“  produced  during  a Hoi  1 ini?  test.  As  a result,  another  method  was  developed  to  remove  the  gas 
from  the  liquid.  It  was  found  that  the  carbon  tetrachloride  could  be  degassed  satisfactorily  by 
evaporating  the  liquid  and  then  recondensing  the  vapor  while  drawing  a high  vacuum  to  remove  the 
noncos-dessablc  saseSi.  This  tedious  process  was  carried  out  using  the  system  shown  in  Figure  34. 
All  containers  and  iineB  were  fabricated  of  Byrcx  glass.  The  procedure  for  negusaiuit  was  as 
follows: 


The  carbon  tetrachloride  was  poured  into  container  A,  and  container  B was 
surrounded  by  a liquid  nitrogen  bath.  Heat  was  supplied  to  container  A to  cause 
the  liquid  to  evaporate.  The  vapor  flowed  over  to  container  B where  it  condensed 
and  solidified.  The  vacuum  pump  maintained  the  system  at  about  2 inches  Hg  ahs 
and  drew  off  the  gas  that  would  not  condense.  After  all  of  the  liquid  was  evaporated 
from  Container  A,  the  liquid  nitrogen  bath  was  removed  from  B and  placed  under  A. 
The  liquid  was  then  evaporated  from  B back  into  A after  suitably  adjusting  the 
control  valves.  The  procedure  had  to  be  repealed  about  four  times  before  the 
desired  amount  of  gas  was  removed  from  the  liquid.  Container  B was  then  filled 
with  the  degassed  liquid  and  valves  1 and  2 closed  and  valves  3 and  4 opened. 
Container  B was  heated  to  increase  the  pressure  and  the  liquid  was  then  drawn 
into  the  test  3cciion.  Figure  34  shown  ...  ■ liquid  level  when  the  system  is  ready 
for  the  start  of  the  test.  Valve  4 was  closed,  and  the  liquid  in  container  B was 
heated  to  give  the  desired  pressure  while  the  liquid  in  the  test  section  was  cooled. 

In  this  manner,  it  was  not  necessary  to  expose  the  degassed  liquid  to  the  atmos- 
phere or  other  contaminating  gas  ir.  order  !o  ohtain  the  required  1 atm  abs  pressure. 

The  test  section  was  similar  to  the  v.atev  test  tank  with  the  exception  that  the  entrance  and 
exit  lines  were  both  glass  welded  to  the  top  of  the  cylinder  body  (Cf.  Fig.  3b/. 

To  maintain  the  carbon  tetrachloride  at  a constant  temperature  during  the  high-subcooling 
tests,  it  was  necessary  to  surround  the  tank  with  a coolant.  Figure  36  is  a schematic  drawing  of 
the  test  tank  surrounded  by  ethyl  alcohol  that  was  used  as  the  cooling  liquid,  A clear  frost-free 
view  of  the  boiling  was  possible  with  use  of  purging  tubes  bled  witl  dry  nitrogen  gas. 


4.  Wall~Temperature  Measurement.  It  was  desirable  to  measure  the  wall  temperature 
required  to  initiate  boiling  in  the  carbon  tetrachloride;  the  excess  wall  temperature  required  for 
water  was  known  (LI.  e.g.  Ref.  15).  The  high  ui;  voltage  gradient  in  the  heeling  strip  "dried  in 
the  usual  difficulty  of  obtainiug  accurate  measurement  of  wall  temper. .dr:-  by  the  direct  attachment 
of  a thermocouple  to  the  wall.  The  most  satisfactory  arrangement  that  was  obtained  for  determining 
the  wall  temperature  is  shown  in  figure  37.  The  hcaliug  ^ trip  was  sealed  Into  tbe  Plexiglass  holder 
with  acetone.  The  thermocouple  assembly  was  screwed  into  the  holder  unlii  electrical  contact  wus 
obtained  between  the  strip  and  the  thermocouple.  The  screw  was  then  rotated  so  that  the  thermo- 
couple was  withdrawn  about  0.0005  inch  from  the  heating  strip.  The  Plexiglass  holder  was  used 
only  during  the  tests  conducted  to  measure  the  wall  temperatures. 


B.  Exp*,r inwntol  Data 

1.  Degassed  Water.  It  was  found  that  the  peak  heal,  flux  data  fer  degassed  distilled 
water  containing  less  than  0.2  cc  air  per  liter  of  water  were  reproducible  within  about  ±10  per 
cent.  In  view  of  the  inherent  statistical  nature  of  the  burn-out  point,  ii  is  not  surprising  that  the 
data  were  not  reproducible  within  closer  limits  (Cf.  p.  M Part  l).  T’nc  values  of  n«*»k  heat  flux 
plotted  in  Figure  38  varied  from  5.8  Btu/'aq  in.-aec  at  80°F  to  0.3  Btu/sq  in.-sec  at  the  saturation 
temperature. 

Employing  the  values  of  wall  temperature  plotted  in  Figure  39,  the  heat  transfer  coefficient 
was  determined.  It  is  seen  from  Figure  40  that  the  hetH  transfer  coefficient  decreased  with  an 
increase  in  liquid  temperature  although  at  a slower  rate  than  the  peak  heat  flux. 
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Tt  ahmild  he  recalled  that  the  liquid  temperature  was  measured  at  a position  about  M such 
directly  above  the  heating  strip. 

A second  thennometer  located  below  and  to  the  side  of  the  strip  measured  temperatures 
about  25”F  lower  at  a heat  flux  of  5 Rtn/si}  iu.-ser  and  10°F  at  a heat  flux  of  1 in. -see. 

The  thermometer  located  /s  iuch  uuove  the  s tr: d has  been  used  as  the  reference  since  its  readine 
probably  was  more  representative  as  an  effective  liquid  temperature. 


2.  Aerated  Water.  Aerating  the  water  so  that  it  contained  more  than  30cc  air  per  liter 
of  water  resulted  in  the  formation  of  hubbies  which  contained  appreciable  quantities  of  noncon- 
dcnsable  gancs.  The  effect  of  the  gas  was  to  lower  the  values  of  the  peak  heat  flux  about  25  per 
cent  below  the  degussed  water  data  (Cf.  Fig.  38).  At  liquid  temperuturea  exceeding  110°F,  a 
large  stagnant  gas  bubble  having  a diameter  about  10  times  as  large  as  the  normal  nucleate  bubble 
would  form  on  the  heater  strip.  This  large  gas  bubble  limited  the  transfer  of  heat  from  the  strip  into 
the  water  and  resulted  in  causing  the  wall  to  melt  at  lower  values  of  heat  flux.  The  peak  heat  flux 
values  that  were  obtained  with  large  gas  bubbles  on  tho  Kent  transfer  surface  also  are  shown  in 
Figure  38.  The  large  gas  bubble  could  be  obtained  occasionally  at  liquid  temperatures  below  11(PF. 

3.  Aerosol-Water  Solution,  if  whs  desir'd  to  investigate  the  effect  cf  a change  ie 

surface  tension  on  the  value  of  peak  heat  flux  without  any  appreciable  change  in  the  other  physical 
properties  of  the  water.  A 1 per  cent  solution  of  Aerosol  (dihexyl  sodium  sulfosuccinate)  and  water 
was  used  since  its  surface  tension  was  only  34  dynes/cm  compared  to  72  dynes/cm  for  distilled 
water.  The  interesting  result  was  that  the  peak  heat  flux  decreased  to  about  half  the  value  for 
degassed  distilled  water  at  80^  (Cf.  Fig.  38). 

The  wail  temperature  for  the  Aerosol-water  f, elution  was  about  15° F lower  than  the  values 
obtained  in  distilled  water.  (Cf.  Fig.  39). 


After  degassing  the  solution  in  the  same  manner  as  the  distilled  water,  the  gas  contest  was 
about  0.7  cc  air  per  liter  of  liquid. 


It  is  interesting  to  notice  that  the  degassed  water,  aerated  water,  and  Aerosol-water  solution 
curves  all  approach  a common  value  of  0.3  Btu/sq  in.-sec  at  the  saturation  temperature  (Cf.  Fig.  38). 


4.  Degasar-.J  Carbor.  Tetrachloride,  The  peak  heat  flux  data  for  the  degassed  carbon 
tetrachloride  also  decreased  with  an  increase  in  liquid  temperature  as  shown  in  Figure  41.  The  data 
easily  were  reproducible  within  +3  per  cent.  > 

Measurements  of  the  wall  temperature  indicated  that  the  «*ce*«  in  wall  temperature  over  the 
saturation  temperature  at  the  peak  heat  flux  was  about  15°  F less  than  for  degassed  water  (Cf.  Fig. 

A<i). 


5.  Aerated  Carbon  Tetrachloride.  Aeration  of  the  carbon  tetrachloride  resulted  in  larger 
hubbies  which  contained  appreciable  quantities  of  a noncondensable  gas.  The  value  of  peak  heat 
flux  was  about  bait  as  large  as  the  degassed  carbon  tetrachloride  at  20PF  (Cf.  Fig.  41).  Unlike  the 
aerated  water  only  one  size  of  gas  bubble  could  be  observed  at  any  given  liquid  temperature. 

It  is  noted  from  Figure  41  that  the  degassed  and  aerated  carbon  tetrachloride  curves  approach 
a common  value  of  0.13  Rtu/sq  in.-sec  at  the  saturation  temperature. 

A qualitative  explanation  of  the  relative  position:  of  the  live  peak  heat  flux  curves  will  be 
given  in  Section  V ufter  discussing  the  proposed  mechanism  of  boiling  heat  transfer. 


III.  MEASUREMENT  OF  BUBBLE  GROWTH  AND  COLLAPSE  RATES 

A.  Description  of  Equipment 

1.  water  Tesis.  Tests  with  water  were  carried  out  using  the  Pyres  g'»»«  container  and 
heating  strip  described  in  the  preceding  section  on  peak  heat  flux  measurements.  A General  Radio, 
35-ram,  shutterless  camera  was  used  with  a Kerr-cell,  electro-optical  shutter  (Cf.  Ref.  22).  The 
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ac rr»ce:;  was  operated  iroro  a 14,1/wvoil  power  supply  lor  a maximum  periou  t«i  u.uj  r>cc.  me 
exposure  time  od  each  frame  was  about  2 microseconds.  A relay  system  allowed  the  camera  to 
accelerate  to  the  required  15,000  frames/sec  in  about  1.2  sec  before  actuating  the  Kerr-cell  and 

firing  the  flash  bulb.  WUh  this  arrangement  about  5 feet  of  the  IGG  feet  of  film  were  exposed. 

The  moderate  resolution  of  the  high-speed  camera  prevented  inspection  of  distances 
smaller  than  0.001  inch. 

2,  Carbon  T ■cirac.hlnride.  Tests.,  The  much  lower  growth  and  collapse  rates  of  the 

carbon  tetrachloride  bubbles  were  determined  by  photographing  the  bubble  at  3000  frames on 
ifi-mm  Eastman  Kodak  Super  XX  motion-picture  film  using  a Fastax  camera. 

The  same  degassing  technique  and  equipment  were  used  during  this  phase  of  the  study 
as  were  described  in  the  preceding  section  on  peak  heat  flux. 


B.  Photographic  Results 

1.  Degassed  Water.  It  may  be  concluded  from  Figure  42  that  the  nucleate  bubble 
formed  somewhere  in  the  liquid  within  0.QC1  inch  of  the  heating  strip.  The  bubble  grew  to  a 
maximum  size  and,  for  liquid  temperatures  below  170°F,  collapsed  while  remaining  withinO.OOl 
inch  of  the  wall.  7hen  the  water  temperature  was  between  approximately  170°F  and  saturation, 
the  bubbles  detached  from  both  t!>c  upper  ana  lower  surfaces  of  the  heating  stnp  and  floated  into 
the  bulk  of  the  liquid  before  collapsing  (Cf,  Fig.  43).  For  siu^licity  of  reference,  these  bubbles 
which  condense  (within  the  visible  0.001  in.)  will  be  denoted  as  vapor  bubbles. 

It  was  observed  «t  all  liquid  temperatures  (Cf.  Fig.  44)  that  the  base  of  some  of  the 
vapor  bubbles  separated  fram  the  heat  transfer  surface  just  prior  to  the  collapse  phase  of  the 
bubble  lifetime  cycle.  Since  this  type  of  bubbie  could  b-  seen  on  the  lower  as  well  as  the  top 
surface  of  the  heating  strip,  the  motiou  of  the  bubble  away  from  the  strip  could  not  be  attributed 

x.  xL.  I ..  . r 
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The  effect  of  heat  flux  on  the  bubble  size  and  population  can  be  seen  from  Figure  45, 
which  shows  several  frames  from  each  of  the  four  motion  pictures  taken  of  the  boiling  at  31,  34, 

63,  and  74  per  cent  of  the  peak  heat  flux.  Five  bubbles  were  measured  from  each  cf  these  films 
and  the  data  plotted  in  Figures  46  through  49.  The  average  values  of  maximum  radius  and  lifetime 
are  plotted  iu  Figure  50.  The  average  of  the  bubble  growth  and  collapse  rates,  hereafter  termed 
bubble  velocity , was  obtained  by  dividing  the  maximum  bubbie  radius  by  half  the  lifetime.  It  is 
seen  that  the  three  factors  are  relatively  insensitive  to  the  magnitude  of  the  heat  flux. 

In  the  next  series  of  testa  the  heat  flux  was  set  at  50  per  cent  of  the  peak  heat  flux  for 
g5r  different  bulk  temperatures  (62,  79,  112,  132,  177° F),  end  met;es 

bubbles.  Several  frames  from  three  of  these  films  are  reproduced  in  Figure  51.  The  data  taken 
from  five  bubbles  at  each  liquid  temperature  are  plotted  in  Figures  52  through  57.  The  average  values 
of  bubble  maximum  radius,  lifetime,  and  bubble  velocity  are  plotted  in  Figure  58.  'The  maximum 
bubble  radius  increased  from  .0135  inch  at  62° F to  .022  inch  at  177° F.  The  lifetime  increased  faster 
then  the  radius  and  resulted  in  a decreese  in  bubble  velocity  from  4.1  ft/sec  at  62°F  to  3.3  ft/sec  at 
177°  F. 


The  large  deviation  from  the  average  values  of  maximum  radios  and  lifetime  for  the  bubbles 
should  be  noted  in  Figures  46  through  49  and  Figures  52  through  57.  As  a result  of  these  large  vari 
ntimig,  the  sycracc  values  clotted  in  Ki cures  50  and  58  nhroild  r..-o  i — -l .t.__ 
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about  25  per  cent  from  the  true  values. 

A qualitative  discussion  of  the  effect  of  liquid  temperature  and  heat  flux  on  the  bubble 
dynamics  will  be  given  in  Section  IV. 

2.  Aerated  Distilled  Water.  The  motion  pictures  of  the  nucleate  boiling  of  aerated  water 
indicated  that  the  bubbles  did  not  condense  while  remaining  against  the  wall  as  was  usually  the 
case  for  the  vapor  bubbles  at  liquid  temperatures  below  170° F. 
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Two  types  of  eaS'vapor  bubbles  were  observed  in  the  serated  water.  One  type  grew  to  about 
the  same  size  so  the  pure  vapor  bubble  end  then  collapsed  leaving  a **!?•“>>  residual  ms  huhhle:  the 
other  type  of  gna-vspor  bubble  grew  to  about  ten  times  the  maximum  size  of  the  vapor  bubble  (*«d 
remained  substantially  stagnant  on  the  beat  transfer  surface.  The  latter,  stagnant  gas-vapor  bubble, 
OvCiiiTCii  most  ircquCQiiv  in  water  eaitL  was  Bi-pristst ur«?ed  with  ^is  and  at  iemperatuifc  above 
iiu°F.  The  collapsing  small  gas-vapor  bubble  secured  most  frequently  in  cooler  water  tfcst  was 
under-saturated  with  gas. 


The  stagnant  gas-vapor  bubble  grew  to  about  the  name  size  ss  the  vapor  bubble  at  approxi- 
mately the  same  rate.  After  this  radius  was  reached,  the  bubble  continued  to  ffow  at  a much  lower 
rate  to  sizeB  up  to  ten  times  the  vapor  hubhle  maximum  radius  (Cf.  Fig.  59).  At  low  heat  flux  several 
of  these  large  bubbles  could  be  seen  cn  the  .trip.  As  the  heat  flux  and  consequently  the  bubble 
population  was  increased,  the  interference  of  neighboring  bubbles  occasionally  would  cause  these 
ras-vnnor  bubbler,  to  detach  from  the  strip  after  reaching  about  the  same  size  ss  the  vapor  bubble. 
However,  usually  the  stagnant  gas-vapor  bubbles  would  cause  burnout  of  the  heater  strip  before 
sufficient  turbulence  could  be  created  in  the  liquid  boundary  to  detach  the  bubble. 


. The  large  gas-vapor  bubble  wouid  decrease  less  than  ten  per  cent  in  radius  when  the  electrical 
power  was  no  longer  supplied  to  the  heating  strip.  This  large  residual  volume  of  gas  indicated  that  the 
bubble  must  have  been  filled  mainly  with  a noncondensable  gas  i ether  than  water  vapor.  The  rate  at 
which  the  gas  in  these  stagnant  bubbles  rediffused  back  iafe  the  liquid  after  the  electric  power  was 
shut  off  was  measured  to  be  less  than  10-5  times  the  rate  at  which  the  bubble  formed. 


The  influence  of  the  liquid  lempcrntsffe  end  hest  flux  nr*  the  dynamics  of  the  short  life  bsb- 
vapor  bubble  was  similar  to  the  effects  on  the  vapor  bubble.  However,  the  large  stagnant  bubble  was 
relatively  unaffected  by  liquid  temperature. 

The  maximum  bubble  radius  and  bubble  velocity  could  not  be  measured  for  heat  fluxes  exceeding 
approximately  25  per  cent  of  the  peak  hest  flux.  At  high  values  of  hest  flux  the  stagnant  bubbles  re- 
maining »n  the  heating  strip  and  the  non  con  den  sing  gas  bubbles  that  detached  from  the  strip  obscured 
the  field  of  view.  The  meager  data  that  were  obtained  indicated  that  the  bubbles  grew  in  the  aerated 
water  at  about  the  same  rate  as  in  the  degassed  water. 


Two  additional  phenomena  were  noticed  in  the  motion  pictures  of  the  aerated  water.  The  first 
was  that  occasionally  the  small  gn« -vapor  bubbles  were  forced  an  appreciable  distance  away  from  the 
heating  surface.  *?hca  these  bubbles  were  ejected  from  the  lower  surface  of  the  strip,  they  traveled 
about  l'A  inches  before  changing  direction  to  hum  • ,.a*J  with  th.-  buoyant  fxce  (Cf.  Fig.  CO).  The 
second  fact  that  could  be  see'j  in  the  motion  pieties,  .ad  a stream  of  warm  liquid  flowing  downward 
against  buoyant  force  from  the  poie  of  the  large  gas  bubbles  cn  the  underside  of  the  strip.  This  warm 
stream  changed  direction  shout  1 inch  helow  the  strip  and,  ««  would  he  expected,  floated  upward. 
G^upliitc  powder  was  placed  on  the  surface  of  these  streaming  bubbles  to  trace  the  liquid  currents.  A 
close  photographic  inspection  indicated  that  the  liquid  film  forming  the  bubble  boundary  was  rotsting 
as  indicated  in  Figure  61. 


These  detaching  and  streaming  bubbles  will  Le  discussed  further  in  Section  IV. 

3.  Aerosol-Water  Solution.  The  nucleate  bubbles  that  were  formed  in  the  degassed  solution 
of  1 per  cent  Asrcnv,-'  and  water  were  generated  within  0.001  inch  of  the  «t?ip  nnrj  grew  obont  the 
same  size  as  the  vapor  bubble.This  type  of  bubble  differed  from  the  vapor  bubble  in  three  main  factors; 
(1)  the  bubble  usually  detached  from  the  atrip  before  collspsing  completely  (Cf.  Fig.  62);  (2)  the  rate 
of  growth  was  considerably  lower  U.uu  for  the  vapor  babble  (Cf.  Fig.  63);  (3)  for  low  population,  about 
50  per  cent  of  the  bubbles  grew  to  the  maximum  radius  and  then  remained  stagnant  on  the  strip. 


The  bubble  maximum  radius  and  lifetime  increased  with  an  increase  in  liquid  temperature  in  s 
similar  manucr  to  the  degassed  water.  The  bubble  velocity  decreased  Slightly  with  an  increase  in 
liquid  temperature  as  seen  in  Figure  64.  The  effect  of  increasing  the  heat  flux  on  the  bubble  dynamics 
was  to  dislodge  the  stagnant  bubbles  from  the  wall  and  thus  increase  the  average  bubble  velocity 
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(Cf.  Fig.  65). 


4.  Degassed  Carbon  Tetrachloride.  The  nucleate  bubbles  formed  in  the  degassed  carbon 
tetrachloride  were  very  similar  to  the  Aerosol-water  bubble.  The  bubbles  formed  on.  or  very  near, 

•he  strip;  grew  considerably  slower  than  the  watef-vspcr  bubble,  and  detached  from  the  wall  before 
condensing.  No  residua'  gas  bubbles  couio  be  sees.  Occasionally  a bubble  was  forced  away  iron- 
the  bottom  of  the  strip.  At  low  heat  flux  shout  50  per  cent  of  the  bubbles  grew  to  a maximum  size 
and  then  remained  stagnant  on  the  heating  atrip  as  was  observed  also  for  the  water-Aerosol  solution. 

Another  effect  couid  be  seen  in  the  shadowgraph  series  cf  motion  pictures.  Figure  66  shows 
that,  as  the  carbon  tetrachloride  bubble  grew,  a thermal  layer  waa  stretched  over  the  bubble  surface 
and  slowly  approached  the  wall. 


The  influence  of  liquid  temperature  on  bubble  dynamics  can  be  seen  in  Figure  67.  The  daiu 
takes  from  this  «**!*■  are  plotted  in  Figures  68  t'  rough  71.  The  bubble  maximum  radius,  lifetime, 
and  velocity  are  plotted  in  Figure  72.  It  is  noted  that  the  bobble  radius  asd  lifetime  increases  while 
the  bubble  velocity  decreases  for  higher  liquid  temperatures.  This  effect  is  similar  to  the  water- 
vapor  bubble. 


Figure  73  summarizes  the  effect  of  heat  flux  on  bubble  dynamics.  The  large  change  in  bubble 
velocity  between  25  end  50  per  cent  of  peak  heat  flux  was  not  observed  lot  the  nucleate  bubbles  in 
water  but  was  similar  to  the  bubbles  forming  in  the  Aerosol-water  solution.  The  photographs  in  Figraen 
$7  and  75  show  the  effects  of  liquid  temperature  and  heat  flux  on  the  size  and  population  of  the  bubbles. 


5.  Aerated  C.nrbon  Tetrachloride.  The  nueleate  bubbles  that  form.;d  in  aerated  carbon 
tetrachloride  were  aimila*  to  the  stagnant  gis-vapor  water  bubbles.  The  bubbles  did  not  collapse 
unless  detached  from  the  strip  by  neighboring  bubbles.  In  this  case  a small  residual  gas  bubble  cuuld 
be  seen  for  several  seconds  after  the  bubble  detached  from  the  heat-transfer  surface. 


The  same  difficulty  waa  encountered  in  measuring  the  maximum  bubble  radius  and  bubble 
velocity  as  for  the  aerated  water.  The  small  quantity  of  data  which  could  be  obtained  indicated  the 
bubble  grew  slightly  slower  in  the  aerated  carbon  tetrachloride  than  in  the  degassed  carbon  tetrachlo- 
ride. 


IV.  PROPOSED  MECHANISM  OF  NUCLEATE-BUBBLE  GROWTH  AND  COLLAPSE 

A Introduction 

The  photographic  study  of  water  and  carbon  tetrachloride  hae  resulted  in  supplying  sufficient 
information  to  suggest  a mechanism  of  the  dynamics  of  nucleate  bubbles  that  is  presented  in  this 
section.  Some  of  the  mechanics  of  the  bull  bit  foists  tints,  grc,‘*h.  and  mltanaa  are  conclusively  shown 
in  the  motion  pictures;  other  factors  are  dedncible  from  the  photographs;  and  some  of  the  remaining 
details  have  not  been  clarified  by  this  study.  It  will  be  seen  that  it  is  possible  to  separate  the  bubble 
life  cycle  into  several  individual  processes. 

There  are  four  basically  different  types  of  nucleate-boiling  bubbles  which  will  be  discussed. 
These  types  occur  respectively  in 

(1)  Degassed-subcooled  water 

(2)  Degassed  bulk-boiling  liquid 

(S)  Aerated  liquid 

(4)  Degassed-subcooled  carbon  tetrachloride 

B.  Byhbls  Formation 

The  first  problem  encountered  in  understanding  the  Lasic  mechanism  of  bubble  dynamics  is 
the  initiation  of  a bubble.  The  study  of  surface  phenomena  i.»  detr-iune  the  role  of  nuclei  in  the  for 
mation  cf  bubbles  is  beyond  the  scope  of  the  present  investigation.  However,  some  of  the  observations 
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on  nucleution  which  were  obtained  incidental  to  this  investigation  together  with  seme  of  the  we!!* 
ksowc  *«?»«  eemirdins?  the  conditions  necessary  to  rupture  a liquid  will  be  discussed  in  this  Section. 

It  is  believed  that  a nucleus  which  contains  gas  or  vapor,  or  a combination  of  both  must  be 
present  in  the  liquid  in  order  to  initiate  the  bubble,  and  that  the  nucleus  is  held  in  a stable  equilib- 
rium position  by  a solid  particle..  To  support  ibis  conclusion  two  conditions  may  be  investigated: 

(1)  the  formation  of  cavitlee  in  a liquid  by  the  random  motion  of  the  molecules  without  the  aid  of 
solid  particles  ana  (2)  stable  gas-vapor  nuclei. 

The  size  of  a spherical  cavity  that  is  needed  to  form  a bubble  may  be  estimated  by  the 
weii-known  relation 

Pi-Po-f^  (?) 


where  (pj  -po)  is  the  pressure  difference  within  the  cavity  that  is  needed  to  overcome  the  surface 
tension  for  s nucleus  radios  of  r.  For  the  esse  of  cavities  forming  within  a liquid,  by  the  random 
molecular  motion,  the  nuclei  would  contain  only  vapor.  Equation  (7)  then  becomes 

Pv  ~ Po  “ (8) 

Equation  (8)  is  plotted  in  Figure  76.  Thia  curve  repreaenta  the  radius  of  the  nucleus  that  is  required 
in  order  tc  have  the  surface  tension  force  of  the  liquid  film  surrounding  a nudeur  in  equilibrium  with 
the  excess  pressure.  For  a given  excess  pressure,  if  the  nucleus  radius  is  lower  than  the  equilibrium 
radius,  the  capillary  pressure  will  exceed  the  ezeers  pressure,  and  the  nucleus  will  collapse  com* 
pletely.  If  the  nucleus  radius  is  larger  than  the  equilibrium  radius,  the  excess  pressure  will  exceed 
the  capillary  pressure  and  the  nucleus  will  grow  into  an  infinite  radius  bubble.  Thus  for  a given 
liquid  temperature  and  pressure  (i.e.,  a given  pv  — p0)  only  one  size  of  vapor  nucleus  may  be  in  equilib- 
rium, and  this  equilibrium  would  be  unstable. 

With  the  aid  of  the  ther!uo^yne~»ic  relation  of  Clapeyron-Clauaiua,  neglecting  the  specific 
volume  of  the  liquid 

dp  A 
dt  " T Pv 

Equation  (8)  then  becomes  for  a perfect  gas,  approximately 


T ~T 

- V 1sst 


..illsatjv  ln/,+2r  \ 


where  Tv  is  the  saturation  temperature  of  the  vapor  within  the  nucleus  at  pressure  py  and  T.«,  is 
the  saturation  temperature  at  pressure  p„.  It  is  known  irom  experimeu  te  that  a heutrtrenafer  surface 
temperature  of  30PF  above  the  saturation  temperature  of  the  liquid  will  esuse  bubble  formation  in 
water.  Substitution  of  this  excess  temperature  into  Equation  (9)  results  in  a required  nucleus  radius 
of  approximately  1ft-4  inch  in  order  to  have  the  pressure  and  surface  tension  forces  in  equilibrium. 
This  value  is  about  J.G4  times  larger  than  the  maximum  cavity  size  that  would  be  expected  to  occur 
by  molecular  fluctuations  with  any  reasonable  frequency  (Cf.  Ref.  23).  As  a consequence,  any  cavity 
that  does  form  in  the  liquid  would  be  smaller  than  the  equilibrium  radius,  and  thus  would  be  col- 
lapsed by  the  surface  tension  force.  The  number  of  nuclei  having  a diameter  exceeding  10"4  in. 
which  arise  from  random  molecular  fluctuations  can  be  estimated  according  to  Volmer  (Cf.  Ref.  27) 
to  be  much  less  than  oue  per  cubic  inch  per  hour. 

Frora  the  above  discussion,  it  is  reasonable  to  eliminate  free  vapur  uuclci  which  arise 
fr  jm  molecular  fluctuations  as  likely  nvcleation  cavities.  We  now  .my  tons  Her  the  possibility  of 
a gas  nucleus  acting  as  a bubble  generation  cavity, 

Seine  brief  mention  of  the  way  in  which  a gas  may  be  held  in  a liquid  will  be  of  interest. 
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A stagnant  liquid  exposed  to  a gas  will  absorb  the  gas  in  the  form  of  single  molecules  dispersed 
lluougfa  the  liquid.  Henry  s Law  states  that  for  a gas  mixture  such  e.s  air,  the  number  of  single 
molecules  of  any  of  the  component  gases  in  solution  with  the  water  will  vary  as  the  solubility  and 
partial  pressure  of  the  component  gas.  When  the  liquid  is  agitated,  for  example,  by  pouring  from 
ooe  container  into  another,  submscroscopic  bubbles  of  the  surrounding  gas  may  be  trapped  ic  the 
liquid  to  form  a colloidal  solution.  If  the  liquid  is  not  completely  saturated  (i.e.  filled  with  the 
maximum  number  of  single  gas  molecules),  these  gas  nuclei  will  eventually  dissolve  into  the 
liquid  to  epproach  a saturation  condition.  This  rate  of  saturating  has  been  estimated  by  Epstein 
and  Plesset  (Cf.  Ref.  24)  to  be  of  the  order  of  several  seconds  in  waier  which  is  undematurated 
with  the  gas.  However,  the  gas  nuclei  may  be  stabilized  on  nonwettable  solid  particles,  and  exist 
indefinitely  in  the  liquid.  In  addition  to  the  gas,  these  nuclei  would  contain  vapor  at  the  satura- 
tion pressure  corns)  jading  iu  ihc  liquid  temperature,  Tee  gnn-v«por  cavities  Would  lie  surrounded 
by  the  surface  of  the  solid  and  a convex  surface  of  the  liquid.  For  this  condition,  the  total  pressure 
within  the  cavity  would  he  less  than  the  liquid  pressure.  The  partial  pressure  of  ihc  gas  would  be 
reduced  below  the  pressure  corresponding  to  a saturation  condition  of  the  liquid,  sad  as  a result 
the  gas  would  no  longer  continue  to  diffuse  into  the  liquid. 

To  see  bow  these  gas-vapor  nuclei  that  are  stabilized  on  solid  particles  may  act  as  gener- 
ation points  for  the  Hubbles.  consider  the  effect  of  increasing  the  temperature  of  the  liquid  in 
contact  with  the  nuclei.  The  temperature  of  the  liquid  may  increase  until  the  partial  pressure  of 
the  vapor  within  the  nuclei  increases  to  the  value  required  to  make  the  total  pressure  sufficiently 
great  to  iorce  some  of  the  gas-vapor  mixture  out  of  the  solid.  Assuming  the  gas-vapor  nuclei  are 
spherical  when  forced  out  of  the  solid,  we  raay  rewrite  Equation  (v)  as 

Pv  + Pg-Po  - y- 

ivhsre:  the  sum  of  the  oartial  pressure  of  the  vapor  and  the  gas  (p..  + p_)  is  the  total  pressure  within 

• • - - - » • g = 

the  cavity.  For  a perfect  gas  the  relation  may  be  written 

U.  R T 

Pv-Po-f--^-  (1C) 

r 4/3  n r3 

It  is  seen  that  the  gas  decreases  the  effect  of  surface  tension  and  thus  decreases  the 
critical  excess  pressure. 

Equation  (10)  may  be  rewritten  (with  the  same  approximations  as  before)  in  terms  of  the 
excess  temperature  as 

«T  T / ..\ 

Tv-Ts«t-— ^ — - !n  U + — 1 I (11) 

\ pr  ) 


or  for  small  temperature  differences 


Tv-T*, 


2a„,fRT:„. 


- a - — 
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The  critical  radius  (i.e.  maximum  radius  for  which  a nucleus  containing  a given  quantity  of 
gas  will  be  in  stable  equilibrium)  becomes 
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(15) 


and  the  excess  pressure  thai  would  be  required  to  cause  the  nucleus  to  form  a growing  bubble 
would  be 


(Pv-Pol 


8tto 

W8R8Tv 


The  term  bubble  will  be  used  to  define  a cavity  having  a radius  greater  than  the  critical  radius; 
the  term  nucleus  will  be  used  to  define  a cavity  having  a radios  smaller  than  the  critical  radius. 

Equation  (10)  is  plotted  in  Figure  77.  It  is  seen  thai  unlike  the  pure-vapor  nuclei,  a 
nucleus  smaller  than  the  ermilihrinm  size  (and  helow  the  critical  radius)  will  ."TOW  to  the  equilibrium 
radius  while  a gas-vapor  nucleus  exceeding  the  equilibrium  si zz  (but  below  tLc  critical  radius)  will 
decrease  to  the  equilibrium  radius.  It  also  may  be  noted  from  Figure  77,  that  decreasing  the  mass 
of  gas  within  the  nucleus  increases  the  value  of  tie  excess  pressure  at  the  critical  radius. 

After  several  hubbies  have  been  formed  from  a nucleus  which  is  stabilized  on  a solid 

naitir-io  mn.it  rtf  tk.  m.  M'nnlii  Km  nitroprl  nn t nf  (Km  fiAVttV  Tn  thi«  nn.n  the  Mnfrntino  flltrUtlH 
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would  contain  vapor  only.  When  this  rondition  is  reached  the  vapor  nuclei  which  have  radii  exceeding 
the  critical  value  when  forced  out  of  the  solid  wiii  form  hubbies;  the  nuclei  which  have  radii  smaller 
than  the  critical  will  collapse  after  leaving  the  solid.  Another  possihle  means  hy  which  pure- vapor 
nuclei  may  form  is  hy  the  random  fluctuations  of  the  liquid  molecules  as  discussed  above.  The 
vapor  nuclei  formed  in  this  manner  and  not  stabilized  on  solid  particles  were  shown  to  be  too  email 
to  form  a bubble  at  any  reasonable  liquid  temperature.  However,  if  the  vapor  nuclei  are  stabilized 
on  solid  particles,  it  theoretically  would  be  possible  for  these  nuclei  also  to  act  as  bubble  gener* 

a»Inn  pclsiS.  Thf  f/jarwofinn  of  KnKKU*  from  yanoa  nna  ja  j « a im|  i y ajwaa  » ♦ m m a 

improbable  condition  of  having  tbe  stabilizing  particle  at  the  point  of  cavity  formation. 


It  is  concluded  that  for  a given  value  of  liquid  superheat  the  gas-vapor  nuclei  which 
exceed  rc  after  being  forced  out  of  the  solid  will  form  bubbles.  As  the  heat-transfer  surface  temper- 
ature increases,  the  vapor  pressure  increases  and  consequently  there  will  be  more  active  or  bubhle 
forming  nuclei  (Cf.  Fig.  77).  Tbe  nuclei  which  contain  the  "mallest  “mount  of  gas  will  have  the 

gT^otoot  £aC6oS  pressure  aud  tcmpcruturc* 

In  addition  to  the  nuclei  that  ere  held  is  the  bulk  of  the  liquid  hy  ouhmicrcscopic  foreign 
particles,  there  would  he  some  gas  trapped  within  the  heat  transfer  surface.  The  capacity  of  the 
foreign  particles  on  the  heat  transfer  surface  to  act  as  nucleation  points  is  determined  mainly  by 
the  cavitv  sizes  and  nonwettiug  character. 


The  buhhles  may  form  more  readily  against  the  heat-transfe>  surface  than  in  the  hulk  of 
the  liquid  if  the  surface  free  energy  of  the  liquid  in  contact  with  vapor  is  greater  than  the  energy 
of  the  liquid  in  contact  with  the  wall.  An  additional  factor  favoring  rapture  of  the  liquid  at  the 
wall  rather  than  in  the  bulk  of  the  liquid  is  the  higher  temperature  at  the  wall.  It  is  evident  that 
»be  surface  condition  of  the  heating  strip  should  affect  the  conditions  necessity  for  nucleation. 


Assuming  the  existence  of  gas-vapor  nuclei  stabilized  on  solid  particles,  it  is  possible 
to  explain  the  high  superheat  that  is  obtainable  in  a liquid  that  has  been  subjected  to  a high 

pnaao;ie# . 0/A  may  inyffihgAfr  th*  fff.i’U  of  inrrfAfting  tit#  | pf???!!??  frx  aayasnl  jtijn/La/) 

atmospheres  and  then  carefully  decreasing  the  pressure  to  the  original  one  atmosphere.  From 
Henry's  Law  it  is  known  that  the  solubility  of  the  gas  in  the  liquid  would  increase  for  higher 
pressure*.  Thus,  some  of  the  gas  within  the  cavity  mould  diffuse  into  the  liquid.  When  the  liquid 
pressure  is  decreased  to  the  original  one  atmosphere,  the  liquid  may  stili  hold  the  gas  in  a super- 
saturated condition.  In  this  case,  tbe  amount  of  gas  in  the  nucleus  would  be  less  than  the  original 
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value.  Thus  when  the  gas-vapor  mixture  is  forced  out  of  the  solid  to  form  a nucleus,  the  partial 
pressure  of  the  gas  would  be  lower.  As  a result  the  vapor  pressure  (superheat)  that  is  required  tc 
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This  high  superheat  hue  been  noted  by  several  investigators  (Cf.  e.g.  Ref.  17). 

Assumptions  that  will  be  necessary  for  discussing  the  generation  of  huhbles  have  been 
presented  briefly  in  tb i **  ni .c mnn <hm:  they  may  be  BunjiuariZcd  as  xollows: 


1. 

2. 


3. 


Gas-vapor  nuclei  stabilized  on  solid  particles  are  present  even  though  the 
liquid  is  not  saturated  with  the  gas. 

The  excess  pressure  necessary  to  form  a bubble  from  a nucleus  which  has 
been  forced  out  of  the  stabilizing  particle  varies  with  the  surface  tension, 
the  partial  pressure  of  the  gas  within  the  nucleus  and  the  nucleus  size 
(f'.f.  Fit.  71. 

' — l ■ *- 

The  liquid  may  ruptuire  near  the  heat  transfer  surface  more  easily  than  in 
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The  Dynamics  of  Bubble*  in  Dscasssd  Subcooled  Water 

1.  Subtle  Gravih  and  Collapse  Mechanism.  The  dynamics  of  nucleate  buhbles  forming 
in  degassed  distilled  water  will  be  the  first  type  discussed  because  these  bubhles  have  the 
simplest  itie  cycle. 

Several  frames  from  a motion  picture  of  the  growth  and  collapse  of  s typical  water  vapor 
huhbSe  are  shewn  in  Figure  42.  The  following  enumerations  refer  to  the  sketches  in  Figure  78. 

The  horizontal  heat  transfer  surface  is  immersed  in  degassed  distilled  water. 

(a)  The  wall  temperature  and  the  adjacent  liquid  film  are  above  the  liquid  saturation 
temperature. 

(h)  A gas-vapor  nucleus  has  been  forced  out  of  a solid  particle.  !f  the  excess  temper- 
ature of  the  liquid  surronnding  the  nuclei  is  greater  than  the  value  given  in  Equation  (13),  the 
equilibrium  pressure  force  would  exceed  the  surface  tension  force  and  the  nucleus  will  become 
a hubble. 

The  motion  pictures  of  the  bubble  formation  indicated  that  the  nucleus  was  within  0.001 
in.  of  the  wall.  For  atmospheric  pool  boiling,  assuming  the  bubble  forms  on  the  wall  which  has  a 
temperature  of  250°F  (Cf.  Fig.  39),  the  excess  pressure  would  be  15  psi.  Thus  from  Figure  77,  it 
is  seer  that  any  nucleus  having  a radius  greater  than  3xl0_v’  in.  wiil  grow  to  form  a bubble.  This 
estimate  is  in  agreement  with  Parkin’s  prediction  (Cf.  Ref.  25)  thst  nuclei  having  a radius  of  the 
order  of  10‘  5 in.  would  initiate  cavitating  bubhles.  Every  hubble  that  forms  at  a wall  temperature 
of  250°F  would  grow  from  nuclei  containing  more  than  2.3  xlO-*®  lbs  of  gas.  For  a constant 
buLLle-Iiijuid-f.lm  temperature,  uic  rale  uf  fpOWlu  of  trie  bubbles  will  depend  0*1  the  difference 
between  the  excess  pressure  existing  at  the  critical  radius  and  the  equilibrium  excess  pressure 
(Cf.  Fig.  77)-  However,  since  the  liquid  film  surrounding  the  bubble  will  be  cooled  by  the  latent 
heat  removed  in  forming  the  vapor  and  the  heat  transferred  to  the  hulk  of  the  liquid,  the  growth 
rate  wilt  he  lower  than  the  isothermal  rate  (Cf.  Fig.  77). 

(c)  The  resistance  of  the  wall  to  the  spherical  growth  of  the  bubble  results  in  approxi- 
mately a hemispherical  shape  (Cf.  Fig.  42).  It  is  probable  that  the  bottom  of  the  buhhle  is  sepa- 
rated from  the  wall  by  a thin  layer  of  liquid  since  the  viscous  <u*.d  adhesion  forces  in  the  liquid 
prevent  the  hubble  from  wiping  the  wall  clear  of  water  and  the  evaDoretion  rate  mav  not  be  suf- 


ficiently large  to  evaporate  any  appreciable  thickness  of  liquid  at  the  bubble 
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against  the  thin  layer  of  water  next  to  the  wall  during  the  growth  portion  of  the  lifetime  by  the 
inertia  force  of  the  pool  of  water.  It  was  aot  possible  to  verify  the  existence  of  this  thin  film 
because  of  the  limited  resolution  power  of  the  camera  lens  system.  However,  the  proposed 
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mechanism  of  growth  or  collapse  of  the  bubble  will  be  the  seme  regardless  of  the  existence  or  non- 
existence of  this  him. 

The  layer  of  sujrerheated  water  thei  was  above  the  nucleation  volume  is  pushed  away  from 
the  wall  by  the  bubble  (Cf.  shadowgraph  Fig.  66).  This  water  film  is  titretched  over  the  top  of  the 
bubble  with  the  result  that  the  temperature  gradient  within  the  film  is  increased.  Heat  is  removed 
from  the  displaced  film  by  conduction  and  convection  to  trie  cooler  surrounding  liquid  ^dotted 
arrows)  and  by  evaporation  of  liquid  into  the  bubble  (solid  arrows).  Tbe  probable  shape  of  the 
temperature  profile  in  the  liquid  at  the  top  of  the  bubble  is  shown  in  Figure  79.  It  is  seen  that  the 
heat  would  be  removed  from  the  film  and  transferred  to  the  bulk  of  the  liquid  as  well  as  into  the 
bubble. 

The  liquid  surrounding  the  bubble  haa  been  set  in  motion  by  the  motion  of  the  bubble 
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(d)  At  the  InsUsis  shown  is  this  sketch  the  film  at  the  ion  of  tbe  bubble  has  been 

cooled  to  the  vapor  saturation  temperature,  and  evaporation  ceases  from  this  region.  The  bubble 
continues  to  grow  although  slower  than  previously  because  of  the  lower  rate  of  evaporation  per 
unit  bubble  volume.  The  liquid  film  surrounding  the  bubble  still  loses  heat  by  evaporation  from 
the  entire  bubble  surface  to  the  cooler  bulk  liquid. 

The  excess  pressure  within  a spherical  burble  during  growth  may  be  estimated  by  the 
Bernoulli  relation  which  for  an  inviscid-iccompressible-infinite  liquid  with  zero  surface  tension 
becomes 
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The  capillary  pressure  may  be  neglected  for  a bubble  having  a radius  of  0.005  inch  (pj-p0 
- 2o/r  < 0,2  psi).  Using  a typical  growth  rate  of  10  ft/sec  (Cf.  c.g.  Fig.  49)  for  the  bubble  radius 
from  .001  inch  to  .010  inch  and  neglecting  the  acceleration  term,  Equation  (14)  may  be  used  to 
obtain  an  average  value  for  the  excess  pressure,  as  follows: 
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f1*1"  corresponding  excess  of  tbe  temperature  of  *he  vepor  within  the  babble  is  6°F. 

These  values  indicate  the  order  of  magni‘  rde  of  the  : .:cess  pleasure  and  temperature  that 
would  be  expected  *o  exist  within  tbe  bubble.  An  exact  solution  would  require  corrections  for  iLe 
influence  of  the  wall  and  mass  addition.  For  the  present  discussion  it  is  sufficient  to  note  that 

1C  n.J  ......  n^.an,.  will  in  tK.  initial  cavitv  In  balance  |!>e  surface  tension  force  has  been 

- *v  I r*  ~ * " r 

dissipated  to  2.C  average  pressure  by  the  time  the  bubble  has  reached  .005-inch  diameter.  This 
excess  pressure  of  2 psi  is  the  average  unbalanced  force  which  caused  the  bubble  to  grow.  The 
excess  temperature  of  the  vapor  also  has  decreased  markedly  from  40  to  6°F.  As  suggested  on  page 
22  insufficient  heat  was  conducted  to  the  bubble  liquid  surface  to  maintain  an  isothermal  condition. 

It  is  possible  to  investigate  the  plausibility  of  the  proposed  mechanism  further  by  con- 
sidering the  maximum  evaporation  rate  that  may  be  expected  under  these  > ondinous  arid  comparing 
this  value  with  the  observed  bubble  growth  rate.  Using  an  analysis  similar  to  that  given  in  Reference 
18,  tbe  heat  tbat  was  removed  from  the  liquid  to  form  a spherical  bubble  of  0.005-inch  radius  would  be 
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The  thickness  of  tb?.  liquid  film  from  which  this  heat  may  be  removed  was 

5-  Dt  - v^SxKT4  (SOxlO-6)-  1.1  xKT4  in. 

Hemoving  7.0  xlO~ ' Btu  from  a liquid  film  l.L  xiO~ * inches  thick  would  lower  the  film  temperature 


approximately 
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6°F 


Assuming  no  other  loss  of  heat  and  a straight  line  temperature  gradient,  the  bubble  liquid  surface 
temperature  would  be  [250  - 2 (6)1  - 238°F. 

The  evaporation  iuto  the  bubble  may  be  calculated  from  the  Knutsen  equation  corrected 

W ith  the  W r r-nmnrlnf inn  paa 
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f-ss€^/~  v/t!  * ! = .0d  / 85-*602.2i  .052  * 9.3  ft/sec  (15) 
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In  this  calculation,  a value  of  .04  was  used  for  the  accommodation  coefficient  as  suggested  ir. 
Reference  26. 

In  view  of  the  assumptions  employed  in  this  estimate,  the  close  agreement  between  the 
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obtained  hy  increasing  the  heat  flux,  the  agreement  is  not  as  good. 
Nevertheless,  the  agreement  of  the  rough  calculation  with  experimental  data  even  at  higher  wall 
temperatures  is  sufficiently  close  to  support  the  proposed  mechanism.  For  an  accurate  analysis 
corrections  would  have  to  be  made  for  the  non-spherical  bubble  shape,  the  beat  transfer  to  the 
subcooled  bulk  liquid,  mass  addition,  latest  heat  Uansport,  and  for  an  initial  bubble  film  tempera- 
ture lower  than  the  wall  temperature.  With  these  modifications  the  agreement  between  the  estimated 
bubble  growth  rate  and  the  measured  growth  rate  ja  expected  to  be  close  for  all  wall  temperatures. 

(e)  in  Sketch  e the  water  surrounding  the  top  region  of  the  bubble  is  now  well  below  the 
saturation  temperature  and  some  of  the  vapor  within  fhe  bubble  condenses  at  the  pole.  The  base 
of  *he  bobble  is  still  surrounded  by  superheated  liquid.  As  a result,  there  is  a flow  of  vapor  from 
the  hot  base  of  the  bubble  to  the  cool  pole.  [Toe  transport  of  heat  by  the  flowing  vapor  has  been 


calculated  ; 


s*  (Cfi  Ref,  27)  for  a one  dimcusionol  flow  to  be  of  a lower  order  of  magnitude 


than  the  total  heat  flux  in  nucleate  boiling.  In  addition  this  condition  of  vapor  streaming  thru- 

the  bubble  occurs  for  only  a fraction  of  the  total  time  that  heat  is  removed  from  the  heating  strip. 

Consequently,  the  flow  of  vapor  through  the  bubble  may  be  eliminated  as  a prime  source  of  heat 

— 1 
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The  bubble  shown  iu  Skstcb  e is  drawn  for  the  condition  where  the  evaporation  rate  is 
equal  !"  the  condensation  rate.  Under  this  condition  there  would  be  no  net  flow  of  vapor  into 
the  bubble;  The  liquid  surrounding  the  bubble  still  has  an  appreciable  outward  velocity,  it  is 
this  *uuiu6uiuiT'  that  viao  stored  iu  *hc  liquid  during  the  growth  of  the  bubble  that  tends  to  curry 
the  bubble  away  from  the  heat  transfer  surface  and  into  the  bulk  of  the  liquid  after  the  bubble 
has  reached  a maximum  radius.  The  photograph  in  Figure  44  clearly  shows  the  bubble  detor- 
maticn  that  results  from  resisting  the  drag  of  the  surrounding  liquid.  Since  this  bubble  was  on 
the  bottom  surface  of  the  heating  Rtrip,  buoyancy  cannot  account  for  the  deformation. 
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(0  I"  this  figure  the  inertia  force  of  the  water  has  been  dissipated  and  the  excess  pressure 
within  the  bubble  changed  from  positive  to  negative.  The  bubble  has  been  deformed  and  the  volume 
increseed  by  the  n«»wn»d  flow  of  the  surroundins  liquid.  The  deformation  also  results  in  moving  the 
bubble  base  out  into  the  cooler  liquid.  The  evaporation  from  the  base  of  the  buLLle  decreases  and 
the  condensation  from  the  pole  increases.  !s  this  figure  the  net  condensation  results  in  a negative 
excess  pressure  within  the  bubble  allows  the  water  to  collapse  the  bubble  back  to  the  same  size  as 
in  Sketch  e. 

(g)  However,  cow  the  rates  of  evaporation  and  condensation  are  not  equal  as  they  were  for 
this  radius  during  the  growth  period.  More  heat  has  been  removed  from  the  superheated  film  and 
transferred  to  the  bulk  of  ihe  liquid.  As  a consequence,  condensation  will  occur  over  a larger  area 
and  the  bubble  will  continue  to  collapse  (Cf.  Fig.  42). 

(h)  The  regions  of  evaporation  and  condensation  are  separated  by  the  saturation  temper- 
ature isothermal  indicated  in  Sketch  h.  The  liquid  surrounding  the  top  of  the  bubble  coniiuues  to 
lose  best  is  the  cosier  surrounding  liquid,  usd  the  saturstios  temperature  isotherm  approaches  the 
wall.  The  cooler  surrounding  liquid  now  has  an  appreciable  inward  radial  velocity. 

(i)  If  the  bubble  collapses  faster  than  the  isotherm  can  approach  the  wall,  once  again, 
as  in  the  growth  portion  of  the  life  cycle,  the  rate  of  evaporation  could  equal  the  rate  of  con- 
densation. Even  under  this  extreme  condition  the  bubble  usually  will  not  stop  the  collapse  since 
the  displaced  superheated  liquid  film  continues  to  trfiasfar  heat  is: _ the  suhrooisd  Knit  liquid. 

This  transfer  of  heat  lowers  the  saturation  isotherm  and  again  causes  the  condensation  to  exceed 
the  rvnuurniiuD.  The  inrush  ing  water  also  acts  ms  an  additional  collapsing  force  by  compressing 
the  vapor  within  the  bubble,  thiie  increasing  the  vapor  saturation  temperature  and  resulting  in 
increased  condensation.  This  increase  in  saturation  temperature  would  not  be  expected  to  occur 
until  the  collapse  velocity  of  the  bubble  exceeds  the  condensation  rate  near  the  end  of  the  collapse 
portion  of  the  bubble  lifetime. 

There  is  one  possibility  for  the  bubble  to  stop  the  collapse  and  remain  stagnant.  This 
condition  may  be  reached  for  a suitable  combination  of  'iuuiu  subcooliug,  conductivity,  latest 
heat  and  condensation  rate.  In  this  case,  the  heat  supplied  by  condensation  would  balance  the 
heat  lost  to  the  bulk  liquid.  This  stagnant  bubble  was  never  observed  in  degassed-subcooied 
water. 

(j)  The  inrushing  cool  water  has  lowered  the  saturation  temperature  area  below  the  free 
convection  position,  and  the  bubble  has  collapsed  back  to  the  critical  size. 

(k)  The  hot  wall  is  (.plashed  with  the  cool  liquid  and  the  local  wall  temperature 


decreases. 


may  repeat. 


The  liquid  near  the  wall  has  been  once  again  superheated,  and  the  bubble  cycle 


2.  Bubble  Lifetime  Cycle  on  Excess  Pressure  Curve.  It  is  interesting  to  trace  the  bubble 
growth  and  collapse  cycle  as  described  above  on  a graph  of  excess  pressure  vs.  bubble  radius. 
Figure  77  shows  the  equilibrium  positions  for  a typical  nucleus.  Point  A on  the  graph  is  the 
assumed  starting  equilibrium  position.  As  the  heat  transfer  surface  temperature  is  increased, 
the  temperature  of  the  liquid  film  surrounding  the  nucleus  also  increases.  The  resulting  higher 
vapor  presswe  causes  ihe  nucleus  radius  to  increase  so  that  the  growth  of  the  nucleus  follows 
the  hue  A-B=  when  the  critical  liquid  film  temperature  is  exceeded,  the  nucleus  grows  to  form  a 
bubble.  Since  heat  is  removed  from  the  liquid  film  surrounding  the  bubble,  the  vapor  pressure 
falls  beiow  the  isothermal  line.  The  LiiLh'o  grows  '•••.Uov.'isg  a line  such  as  At  Point  f!. 
sufficient  heat  has  been  removed  from  the  licruid  film  to  lower  the  bubble  path  line  to  the  equi- 
librium position.  The  bubble  would  remain  stagnant  at  this  point  were  it  not  for  the  transfer  of 
heat  from  the  liquid  film  into  the  cooler  hulk  liquid.  Ac  additions!  factor  which  decreases  the 
vapor  pressure  in  ihe  bubble  is  the  decreased  evaporation  as  e result  of  the  movement  of  the 
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bubble  base  away  from  the  hot  wail.  Ae  a result  of  decreasing  the  average  liquid  film  temperature* 
thn  vapor  pressure  within  the  bubble  falls  below  the  suiTounding  liquid  pressure.  The  bubble  then 
collapsed  back  to  radius  A along  a path  such  as  C'D. 

As  soon  as  the  liquid  film  has  received  sufficient  heat  from  the  wall,  the  cycle  may 

3.  influence  of  Liquid  Conditions  on  Bubble  Dynamics.  In  order  to  substantiate  further 
the  proposed  mechanism  of  bubble  dynamics  in  degassed  water,  an  attempt  will  be  made  to  predict 
the  effects  of  liquid  temperature,  wali  temperature,  heat  flux  and  liquid  pressure  on  the  bubble 
size  and  velocity.  These  predicticne  will  then  be  compared  with  the  experimental  data. 

(a).  Effect  of  Liquid  Temperature.  In  the  proposed  mechanism  ri  bubble  growth  it  was 
indicated  that  the  degree  of  superheat  in  the  liquid  film  controlled  the  growth  rate,  Equation  (15). 
Since  the  initial  superheat  is  determined  Ly  the  quantity  of  gas  within  the  nuclei  and  the  magnitude 
of  the  liquid  surface  tension;  Equation  (11).  it  would  be  expected  that  the  inifi"!  jfnwth  rate  would 
remain  substantially  independent  of  iiquid  temperature.  The  average  growth  rate,  however,  would  be 
lowered  slightly  for  a lower  liquid  temperature  because  of  a more  rapid  decrease  of  tie  temperature 
of  the  liquid  film  surrounding  the  bubble.  Th*  collapse  rate  should  be  more  strongly  influenced  by 
the  liquid  temperature*  The  bubble  collapse  rate  would  increase  with  decreasing  liquid  temperature 
because  of  the  iower  bubble  film  temperature  and  consequently  greater  condensation  rate.  These 
effects  are  seen  in  the  experimental  results  shown  in  Figures  52  through  58.  In  these  Figures  the 
life  cycles  of  five  different  babbles  are  shown  for  liquid  temperatures  from  62° F to  177°r.  The 
pressure  was  1 atm  for  all  of  these  cases  and  the  heat  flux  was  at  50  per  cent  of  the  peak  value. 
Under  those  conditions  the  average  bubble  growth-rate  was  within  8 per  cent  of  3.8  ft/sec  for  the 
six  different  temperatures  whereas  the  collapse  rates  were  5.6,  4.0,  3.9,  3.5  and  3.1  ft/sec  at  62, 
78,  125,  135,  177°F,  respectively. 

For  a constant  pressure,  the  bubble  should  reach  a smaller  maximum  radius  as  the  liquid 
bulk  temperature  is  decreased.  THi"  tredieiion  follows  directly  from  the  fact  that  the  growth  rate 
would  be  constant  or  slightly  decreased  and  the  knowledge  that  a lower  bulk  temperature  will 
decrease  the  growth  time  by  removing  the  superheat  from  the  liquid  film  surrounding  the  bubble 
more  rapidly.  Consequently,  the  bubble  radius  at  which  condensation  and  evaporation  are  equal 
should  decrease  for  a lower  liquid  bulk  temperature.  Reference  to  Figure  58  shows  that  the  typical 
maximum  radii  decreased  from  .022  in.  to  .0135  in.  as  the  bulk  temperature  was  decreased  from 
177°F  to  62°F»  The  trend  c£  the  experiment"!  data  agrees  with  the  prediction. 

Since  (1)  the  growth  rate  is  substantially  independent  of  liquid  temperature,  (2)  the 
collapse  rate  increases  for  lower  liquid  temperature,  and  (3)  the  maximum  bubble  radius  decreased 
for  decreased  liquid  temperature,  it  is  obvious  that  the  total  lifetime  should  also  decrease. 
Deferring  again  to  the  results  shown  in  Figure  58,  the  average  lifetimes  of  tho  bubbles  decreased 
from  1100  microsec  at  177° F to  500  microsec  at  62° F. 

As  seen  from  the  foregoing  comparison,  the  growth  rate,  maximum  radius  and  lifetime  as 
determined  from  the  motion  pictures  actually  are  influenced  by  the  bulk  temperature  of  the  liquid 
as  predicted  from  the  proposed  mechanism. 

(u)  Effect  of  Wall  Temperature.  The  temperature  of  the  heat  transfer  surface  may  be 
increased  by  supplying  more  power  to  the  beating  strip.  As  a result  the  adjoining  liquid  film 
temperature  will  increase  and  result  in  activating  a greater  number  of  nuclei.  The  bubbles  will 
create  turbulence  iu  the  liquid  close  to  the  v.' all  ccd  thus  coo!  the  liouid  film  sur*‘r' nnitino  1 1.^ 

» l O O'” 

baring  hubbies.  This  increased  cooling  rate  over  the  value  obtained  at  a lower  bubble  population 
decreases  the  time  needed  to  iower  the  liquid  film  temperature  iu  toe  vaiUc  needed  to  balance  the 
evaporation  and  condensation  rules.  Because  of  the  increased  ceding,  the  collapse  rate  of  the 
bubble  should  increase.  Thus,  an  increase  in  the  heat  flux  (i.c.  wail  temperature)  should  result 
in  lower  maximum  bubble  radius  and  shorter  lifetime. 
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Figures  46  through  49  "hew  buhhle  cycles  at  1 atm*  pressure  and  at  heat  fluxes  of  31,  34,  63 
74  j>sr  cent  of  the  p*nk  heat  flux.  The  correa ponding  maximum  radii  are  0.0195  in.,  0.0195  in., 
0,0185  in.,  and  0.0175  in.;  the  lifetimes  are  900  s.icrosec,  90U  micro  sec,  850  microscc,  end  SQ0 
miciosec,  respectively.  The  decrease  in  size  and  lifetime  that  was  predicted  is  evidenced  by  these 
experimental  data. 

it  should  be  realized  that  the  possible  increase  in  wall  temperature  by  increasing  the  heat 
flux  is  limited  by  the  rapid  increase  in  huhble  population  with  wall  temperature  which  soon  results 
in  hurnout  of  the  heating  strip.  The  effects  discussed  above  were  correspondingly  small.  The  effect 
of  increasing  waii  temperature  while  maintaining  the  bubhle  population  constant  (e.g.  by  increasing 
the  liquid  pressure)  will  be  discussed  in  the  aext  section. 

(c)  Effect  of  Pressure . !t  has  heen  proposed  in  the  preceding  pages  that  a bubble  was 
generated  when  the  vapor  pressure  was  sufficiently  great  io  overcome  the  surface  tension  force. 

The  icu^eraturc  of  the  liquid  film  surrounding  the  hubble  corresponded  to  the  saturation  temperature 
at  the  required  vapor  pressure.  The  growth  rate  was  controlled  by  the  rate  of  evaporation  ists  the 
buhhle.  and  the  evaporation  rate  was  controlled  by  the  liquid  film  temperature.  As  s consequence 
of  the  bubble  generation  bring  controlled  by  the  vap.T  pressure,  a large  change  of  saturation  temper- 
ature with  pressure  would  result  in  & high  hubble  growth  rate.  With  this  fact  in  mind,  we  may  discuss 
the  influence  of  pressure  on  the  bubble  dynamics. 

The  effect  of  increased  liquid  pressure  would  be  to  tend  to  drive  the  gas  in  the  nuclei  into 
solution  with  the  liquid.  However,  the  solubility  may  not  increase  since  the  liquid  film  temperature 
would  increase  with  higher  pressures  and  thus  tend  to  decrease  the  solubility.  Assuming  as  a iirni 
approximation  that  the  mass  of  gas  contained  within  soy  nucleus  does  not  decrease  at  higher 
pressures,  the  critical  excess  pressure  at  a liquid  pressure  of  2000  psi  as  calculated  hy  Equation 
(13)  would  be  5 psi.  The  corresponding  excess  temperature  is  less  than  1°F.  The  excess  pressure 
at  the  high  liquid  pressure  was  uhout  Ms  the  value  at  1 atm  because  the  higher  wall  temperature 
resulted  in  a lower  surface  tension.  The  excess  well  temperature  was  Vi&O  of  the  value  at  1 atm 
because  of  the  small  change  of  saturation  temperature  with  pressure  at  a liquid  pressure  of  2000 
psi. 

From  Figure  81  taken  from  Reference  3,  it  can  be  seen  that  the  measured  excess  temper- 
ature for  degassed  water  decreases  from  60PF  at  1 afen  to  5°F  at  2000  poi.  The  lack  of  solubility 
data  for  water  at  high  pressures  a&d  temperatures  prevents  a rigorous  calculation  from  being  made 
to  predict  the  excess  temperature  variation  with  pressure.  The  ■■•'•motion  of  uo  increase  is 
solubility  at  the  higher  liquid  pressures  results  in  an  estimated  excess  temperature  curve  having 
the  same  trend  but  considerably  below  ibe  experimental  data  (Cf.  Fig.  81)-  If  the  solubility  of  air 
in  water  does  increase  at  higher  pressures,  the  two  curves  would  be  in  better  Agreement. 

Having  estimated  that  the  excess  temperature  would  decrease  with  an  increase  in  liquid 
pressure,  we  may  predict  the  influence  of  pressure  on  the  behavior  of  the  nucleate  buhhles. 

The  lower  excess  temperature  at  higher  pressures  would  decrease  the  evaporation  rate 
and  consequently  decrease  the  huhhle  growth  rate.  This  prediction  is  compatible  with  the  small 
amount  of  data  available  on  the  growth  rate  of  bubbles  at  high  liquid-pressures. 

The  alternate  case  of  boiling  is  r liquid  having  a very  low  pressure  now  may  be  con- 
sidered. The  large  change  of  saturation  temperature  with  pres?!”’**  would  result  in  high  evapo- 
ration rates.  The  bubble  would  be  expected  to  grow  more  rapidly  in  the  liquid  having  a lev? 
pressure.  Ab  an  example,  bubbles  fonuiug  in  water  subjected  to  »n  absolute  pressure  of  0.1 
f*«  woald  be  expected  to  have  a liquid  tiim  excess  temperature  of  150°F  as  calculated  from 
Equation  (11).  'The  high  evaporation  rate  ke  this  excess  temperature  would  maintain  a high 
•Keens  pressure  while  the  hubble  grew  to  a much  larger  size  than  bubbles  forming  in  a liquid 
having  a higher  pressure  and  equal  nubcooling.  The  extremely  high  su-perhest  should  cause  the 
bubhles  to  form  in  an  explosive  manner.  The  well-known  high  bubble-growth  rates  obtained  in 
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vacuum  distillation  apparatus  substantiates  this  prediction. 

The  collapse  rate  of  the  bubble  depends  on  the  liquid  subcooling  and  the  condensation 
rate.  For  constant  subcooling,  tbe  bubble  collapse  rate  should  be  lcwer  at  the  higher  pressures. 

The  change  in  latent  heat  of  evaporation  has  to  be  considered  at  the  very  high  pressures. 
If  the  latent  heat  decreases,  die  heat  which  has  to  be  taken  from  the  liquid  surrounding  the  babble 
tc  fill  the  cavity  with  v«nnr  is  decreased  also.  The  lower  latent  heat  will  tend  to  keep  tbe  temper* 
ature  of  tbe  film  surrounding  tbe  growing  bubble  at  a value  closer  to  the  initial  film  temperature 
and  thus  tend  to  increase  the  evaporation  rate.  The  initial  difference  between  the  film  temperature 
and  the  saturation  temperature,  however,  is  so  strongly  affected  by  pressure,  as  mentioned  above, 
that  the  effect  of  the  latent  heat  should  be  almost  negligible  in  comparison. 

D.  The  Mechanics  of  Bubble  Growth  and  Caiiapse  in  a Degossed>8gik  Boiling  Liquid 

The  formation  and  growth  of  a bubble  in  a degassed*bulk  boiling  liquid  should  be  the  same 
as  in  the  degassed-cubcooled  liquid.  The  oniy  difference  is  that  the  maximum  bubble  size  is  iarger 
for  the  bulk  boiling  as  would  be  expected  from  the  discussion  in  Section  iV-C*3.  Since  no  net 
condensation  can  occur,  the  bubble  will  grow  until  tbe  buoyant  force  or  surrounding  liquid  current 
carries  it  off  of  im>  heat  transfer  surface.  At  high  bubble  population  tbs  neighboring  bubbles  may 
create  sufficient  twbulence  in  tbe  liquid  to  wipe  the  bubble  off  of  the  wall.  Several  frames  of 
hulk  boiling  water  shown  in  Figure  82  support  these  conclusions.  It  should  be  noted  that  bubbles 
lorming  in  a liquid  having  a oumi  i SubcoGirug  wouia  sisu  ctkihu  this  life  cycle.  In  th*.*  case,  the 
condensation  rate  is  small  so  that  tbe  bubble  may  travel  an  appreciable  distance  away  from  the 
wall  before  condensing  (Cf.  Fig.  83). 

The  bubbles  forming  in  a bulk  boiling  liquid  would  be  expected  to  have  a smaller  growth 
rate -at  higher  liquid  presuires  for  the  same  reasons  as  preseuted  in  the  preceding  section. 


E.  The  Mechanism  of  fha  Dynamics  af  Bubbles  in  an  Aerated  Liquid 

The  life  cycle  of  a bubble  growing  in  a liquid  which  contains  a high  concentration  of  gas 
is  considerably  different  from  tbe  vapor  bubble  described  in  Section  IV  C.  This  gas-vapor  type  of 
bubble  occurred  in  tests  of  water  having  more  than  18  cc  air/liter  water  (water  at  212°F  and  1 atm. 
is  saturated  when  13  cc  air  is  dissolved  in  one  liter  of  water). 

As  before,  the  growth  uuu  collapse  will  Kn  ri i sn:sRf d iii  steps.  The  Sketches  of 

Figure  84  correspond  to  the  enumerations  of  ihe  following  paragraphs. 

(a)  Tuc  liquid  next  to  the  wall  is  superheated. 

(b)  A nucleus  of  sufficient  size  is  present  so  that  th?  bubble  may  form.  Because  of  the 
high  concentration  nf  gas  in  the  liquid,  the  nuclei  would  contain  a great*?  gas  pressure...  4«  » 
result  the  wall  temperature  required  to  initiate  bciiing  would  be  expected  to  be  less  than  for  the 
degassed  liquid  (Cf.  Fig.  77). 

(c)  Since  the  solubility  of  most  gases  decreases  with  temperature,  the  liquid  film 
surrounding  the  nucleus  may  be  supersaturated  with  the  dissolved  gas  even  though  the  bulk  of 
the  liquid  is  undersaturated  with  gas.  The  bubble  to  be  described  in  tbe  first  part  of  this  section 
is  grooving  in  a liquid  which  has  been  saturated  with  gas.  As  a consequence,  the  superheated 
boundary  layer  would  be  supersaturated  with  gas.  (A  bubble  growing  in  a liquid  which  is  under* 
saturated  withs-s-  is  discussed  os  Page  29).  The  liquid  also  contains  an  appreciable  quantity 
of  gas  in  the  form  of  sub-microscopic  gas-vapor  nuclei.  It  would  be  expected  that  the  number  of 
these  nuclei  would  be  greater  than  iu  the  degassed  liquid  and  that  there  consequently  would  be 
a greater  number  of  large  nuclei  that  would  form  bubbles  at  a relatively  low  excess  pressure. 

Thus  the  average  bubble  growth  rate  should  be  somewhat  less  than  for  the  vapor  bubble  for  the 
reasons  given  in  Section  SV-B.  As  the  bubble  grows,  the  gas  contained  in  the  liquid  film  sur- 
rounding the  bubble  would  diffuse  into  the  bubble. 
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(d)  Let  us  now  consider  die  inutant  at  which  the  evaporation  and  condensation  rates  are 
This  instant  is  illustrat'd  in  Sketch  d of  F>fl,*,r«  R4-  Unlike  the  pure  vapor  bubble,  there  still 
will  be  a net  mass  addition  into  the  gna-vapor  bubble  because  of  the  diffusion  of  gas.  If  the  rate  at 
which  the  ata  enters  the  bubble  is  sufficiently  gpreat,  the  bubble  will  maintain  its  shape. 

(*)  The  pertiei  pressure  of  the  vapor  in  the  bubble  decreases  since  the  film  has  lost  heat 
to  the  mnier  bulk  liquid.  The  condensation  rate  now  exceeds  the  evaporation  rate.  Since  the  bubble 
base  was  not  dragged  away  from  the  heat  transfer  surface  by  the  surrounding  low-velocity  liquid, 
the  net  condensation  will  be  small.  If  we  assume  a diffusion  rate;  that  is  higher  than  the  net  conden- 
nation  rate,  the  bubble  will  continue  to  grow  slowly.  This  would  be  the  case  for  a high  concentration 
of  gas  iu  iLe  liquid. 

(f)  In  the  next  step  the  saturation  temperature  isotherm  has  reached  an  equilibrium  position 
where  the  heat  lost  by  the  combined  processes  of  conduction  to  the  cooler  liquid  and  evaporation  is 
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hiihliie  conticuen  to  grow  althmidi  at  as  extremely  low  rate  since  gas  continues  to  diffuse  into  the 
bubble.  The  bubble  contains  as  appreciable  amount  ui  gua  a a Las  been  verified  experimentally.  Wbss 
the  electrics!  power  was  shut  off-  the  bubble  decreased  less  than  10  per  cent  in  radius  and  thus  it  was 
indicated  that  the  majority  of  the  bubble  was  tilled  with  a ncs-coudrnsable  gas.  At  first,  it  is  difficult 
to  exploit.  Uit  lirgc  diffusion  rote  because  the  iicrmai  rats  of  diffusion  of  sir  in  water  is  !e®=  ln»n  lO--' 
times  the  rate  at  which  this  bubble  formed.  One  has  to  realize,  however,  that  the  liquid  film  is  super- 
saturated with  gas,  and  that  an  additional  increase  in  rate  over  the  diffusion  rate  in  still  water  is 
possible  since  the  liquid  film  surrounding  the  bubble  was  supplied  continuously  with  nun-degassed 
liquid  due  to  the  induced  currents.  These  currents  were  examined  experimentally.  Close  photographic 
examination  of  the  bubble  indicated  that  the  liquid  film  was  rotating  at  a maximum  velocity  of 
approximately  3 in./sec.  Figure  61  is  a series  of  pictures  taken  of  graphite  traces  on  the  bubble 
surface.  It  is  concluded  that  the  ass-vapor  bubbles  forming  in  a liquid  containing  a high  concentration 
of  dissolved  gases  will  grow  to  a maximum  size  and  remain  substantially  stagnant  on  the  beat  transfer 
surface.  Although  the  initial  growth  rate  of  these  bubbles  is  clone  to  the  growth  rate  for  ihe  vapor 
bubbles,  the  gas-yepor  bubble  will  not  collapse  when  the  net  evaporation  into  the  bubble  ceases. 

The  bubble  "ontinues  to  grow  because  of  diffusion  of  gas  into  the  cavity.  A series  of  photographs 
showing  a bubble  growing  in  water  saturated  with  air  is  given  in  Figure  85. 

Bubbles  forming  in  water  which  is  undersaturated  with  gas  would  be  expected  to  exhibit  a 
i - v - . — "3  — ’ i — . . both  the  degassed  liuuid  and  svs-satvrated  T'h^  nmta^r  nf  ih. 
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within  the  nuclei  for  this  case  should  be  lower  than  in  the  liquid  that  is  saturated  with  gas.  As  a 
result  the  wall  temperature  required  to  iniiiate  boiling  would  be  slightly  hi^ier  than  for  the  satu- 
rated liquid.  One  would  expect  the  diffusion  of  gas  into  the  bubble  to  be  lower  for  the  under- 
saturated  liquid.  In  such  a case  the  condensation  rate  may  exceed  the  gps  dL.usion  rate  and  the 
bubble  would  tend  to  collapse.  A small  residua!  gas  bubble  would  remain  after  the  vapor  has  been 
condensed.  If  the  gas  bubble  is  not  stabilized  on  a nonwettable  solid,  the  gas  would  diffuse  into 
the  liquid  for  the  reasons  presented  iu  Section  IV-B. 

Bubbles  have  been  observed  in  water  containing  less  than  12  cc  air/liter  o:  water.  In  this 
intermediate  aerated  liquid,  the  bubbles  grew  to  approximately  [he  same  size  as  the  vapor  bubble 


ana  cuiiapncu,  leaving  « Small  residual  ^aS  bubble  os  predicted  (Cf.  Fig.  60). 


F.  TK«  Mechanism  of  Bubble  Growth  nnd  Collapse  in  Osqastad-Subcoolsd  Carbon 

Figure  8 6"  shows  a typical  photographic  sequence  of  a bubble  forming  in  degassed-suiw 
cooled  carbon  tetrachloride.  Is  addition  a series  of  Sketches  is  given  in  Figure  86b  to  illustrate 
the  steps  in  the  following  discussion. 

(a)  In  Sketch  a of  Figure  S6  the  temperature  of  the  beat  transfer  surface  is  above  the 
liquid  saturation  temperature. 
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(b)  Assuming  the  same  size  of  nuclei  as  was  required  to  form  a bubble  in  water,  ii  is 
seen  from  Figure  37  that  the  excess  pressure  would  be  about  5 p*w  The  coiresponding  excess 
temperature  is  about  20° F.  This  value  agrees  with  tbs  experiments!  d«tn  presented  in  Figure  39. 

It  is  iuicrcsiiug  to  note  that  as  assumed  nucleus  radius  of  3xi0-:>  in.  results  in  a correct  estimate 
of  wall  temperature  for  both  the  carbon  tetrachloride  and  the  water. 

As  a result  of  the  lower  excess  pressure,  the  initial  growth  rate  of  the  carbon  tetrcchloride 
bubble  would  be  shout  half  of  the  water  bubble  growth  rate.  The  measured  giowth  velocities  for 
carbon  tetrachloride  bubbles  at  atmospheric  pressure  are  shown  in  Figures  68  through  71.  The 
average  bubble  velocities  for  various  liquid  temperatures  are  shown  in  Figure  72.  It  is  seen  by 
comparing  Figures  58  sad  72  that  *br  carhor  tetrachloride-vapor  bubble  velocity  is  about  % as 
pyeat  as  for  the  water-vapor  bubble  -i  the  same  condition. 

As  a result  of  the  slower  bubble  growth  in  carbon  tetrachloride,  the  surrounding  liquid 
would  be  forced  away  from  the  wall  at  a lower  velocity. 

(c)  The  bubble  continues  to  fpow  thereby  stretching  the  superheated  film  of  liquid. 

Heat  is  transferred  from  the  film  to  the  wubcooled  bulk  of  the  li ouid. 

(d)  The  film  surrounding  the  pole  of  the  bu  hole  has  been  Coded  tG  the  ovituTStlOn 
temperature  and  evaporation  ceases  from  this  area. 

(c)  The  saturation  temperature  isotherm  which  separates  the  evaporation  and  con- 
densation areas  has  approached  the  wail  so  that  the  condensation  area  has  increased. 

(0  All  of  the  liquid  that  is  now  evaporated  from  the  lower  area  of  the  bubble  is 
condensed  on  the  upper  area  of  the  bubble.  Since  there  is  no  net  flow  of  vapor  into  the  bubble, 
there  is  no  further  growth  (except  for  the  small  change  in  the  volume  of  the  bubble  that  is 
necessary  to  decrease  the  excess  pressure  to  a zero  value).  The  low  growth  rate  of  the  bubble 
coupled  with  the  lower  thermal  conductivity  has  resulted  in  a maximum  bubble  size  approximately 
one  half  that  for  the  water  vapor  bubble  (compare  Fig.  58  ana  72). 

(g)  In  Section  1V-B,  it  was  seen  that  the  initial  rate  of  growth  of  any  bubble  in  a given 
liquid  depended  on  the  pressure  of  the  gas  within  the  nucleus.  The  low  thermal  conductivity  of 
the  carbon  tetrachloride  is  insufficient  to  transfer  enough  beat  from  die  liquid  film  surrounding 
the  buobie  to  cause  u net  cuudensaiion,  unless  the  base  of  the  bubble  is  moved  away  from  the 
heat  transfer  surface.  As  a result  of  these  two  factors,  it  may  be  concluded  that  die  bubbles 
originating  from  some  nuclei  may  collapse  while  the  bubbles  originating  from  other  nuclei  will 
grow  to  a maximum  size  and  remain  stagnant. 

Figure  86  shows  a bubble  that  has  originated  from  a nucleus  having  s suitable  geometry 
and  containing  a small  quantity  of  gas  so  as  to  result  in  a high  evaporation  rate  and  thus  a re- 
latively high  growi.h  rate.  The  surrounding  carbon  tetrachloride  has  been  forced  away  from  the 
wall  by  the  growing  bubble  with  sufficient  velocity  to  drag  the  bubble  away  from  die  wall  wheu 
the  net  evaporation  ceeees. 

It  will  be  recalled  that  the  w«t*r*vapor  bubble  that  was  discussed  in  Section  IV-C  had 
a similar  tendency  to  detach  from  the  heating  strip.  However,  the  higher  thermal  conductivity 
in  the  water  resulted  in  greater  rate  of  heat  transfer  so  that  the  water  vapor  was  condensed  and 
thus  the  bubble  collapsed  before  it  could  be  carried  away  from  the  wall.  At  very  low  subcooling, 
(he  heat  transfer  from  the  water  film  into  the  bulk  of  the  liquid  was  lower  and  the  bubble  was 
carried  an  appreciable  distance  away  from  the  he-.v-tisnsfer  surface  before  collapsing. 

Figure  88  shows  a bubble  that  has  originated  in  carbon  tetrachloride  from  a nucleus 
with  conditions  which  result  in  a low  evaporation  rate  and  thus  relatively  low  growth  velocity. 

The  surrounding  liquid  has  insufficient  velocity  to  drag  the  bubble  away  from  tbe  hot  wail.  In 
this  case  the  heat  transferred  from  the  liquid  film  to  the  cooler  bulk  liquid  is  balanced  by  the 
heat  supplied  by  condensation  and  conduction  from  the  hot  wall.  As  n result  tbe  bubble  remains 
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»t«gna!!t  on  the  heating  atrip. 

The  audition  of  Aerosol  to  wat*r  drrrensed  the  surface  tension  of  the  solution  to  approxi- 
mately the  same  value  aa  for  the  carbon  tetrachloride.  As  would  be  expected  from  the  above  dis- 
cussion. bubbles  forming  in  Aerosol-water  solution  resembled  the  carbon  tetrachloride  babble  by 
growing  slowly  and  cites  remaining  stagnant  on  the  strip  (Cf.  tig.  62).  It  may  oe  concluded  U.at 
the  thermal  conductivity  of  the  water  ie  also  insufficiently  great  to  result  in  conducting  irons  heat 
away  from  the  bubble  pole  than  is  transferred  ft  the  pole.  As  a resuit  the  bubble  will  act  collapse 
unless  the  base  is  removed  from  the  hot  liquid  adjacent  to  the  heat  transfer  surface. 

The  remainder  oi  this  discussion  will  describe  the  collapse  phase  of  the  life  cycle  for  the 
bubbles  which  are  dragged  awey  from  the  wall. 

(h)  The  bubble  which  detaches  from  the  heat  transfer  surface  floats  into  the  subcooled 
bulk  of  the  liquid  and  is  replaced  at  the  wall  by  cool  liquid. 

(i)  The  vapor  within  the  bubble  condenses  on  the  subccoied  film.  The  latest  heat  from 
thi-;  condensation  is  transferred  to  the  surrounding  subcooled  liquid.  The  rate  of  condensation  is 
lover  than  for  the  water-vapor  bubble  because  uf  the  smaller  subcooling  in  the  bubble  film  as  a 
result  of  the  low  thermal  conductivity.  The  latent  heat  per  unit  volume  of  vapor  also  is  a factor 
controlling  the  bubble  film  temperature.  However,  this  quantity  is  about  equal  for  “stw  and  carbon 
tetrachloride  under  the  existing  conditions. 

(j)  The  velocity  of  the  liquid  nenr  the  wall  has  decreased  and  the  boundary  layer  has  l-oen 

heated. 

(k)  The  boundary  layer  has  been  superheated  so  that  the  bubble  cycle  may  repeat. 

The  carbon  tetrachloride  vapor  bubbles  can  be  seen  several  maximum  bubble  diameters  away 
from  the  beat  transfer  surface,  while  the  water  vapor  bubbles  usually  collapsed  entirely  while  remain- 
ing within  O.GGi  in.  of  the  wall.  The  major  cause  tor  this  different  behavior  is  the  lower  thermal  conductivity 
of  ihe  carbon  tetrachloride  and  the  resultant  longer  collapse  time. 

The  bubbles  forming  in  degassed-subcooled  carbon  tetrachloride  would  be  affected  by  changes  in 
the  bulk  temperature,  wall  temperature,  and  liquid  pressure  in  a similar  manner  to  the  bubbles  forming  in 
degassed-subcooled  water  for  the  same  reasons  as  given  in  Section  IV-C. 


Y.  PROPOSED  MECHANISM  OF  NUCLEATE-BOILING  HEAT  TRANSFER 


A.  introduction 

The  bubble  growth  and  collapse  cycle  discussed  in  the  preceding  Sections  may  be  used  to 
predict  the  variation  of  boiling  heat  transfer  for  various  liquid  cond>‘ions. 

It  was  discussed  in  Section  IV-C-1  that  the  ret*  of  transport  of  latent  heat  through  the 
bubbie  »«»  insignificant  compared  to  the  measured  total  heat  transfer.  Jakob  (Cf.  Ref.  7)  calculated 
the  heat  required  to  fill  nucleate  bubbles  with  vapor  and  concluded  that  this  mechanism  also  could 
not  account  for  the  high  rates  of  boiling  heat  transfer.  As  a result,  jakcb  proposed  that  the  high 
rates  of  heat  transfer  that  were  observed  in  pool  bulk  boiling  were  the  result  of  the  agitation  of  the 
liquid  !!t=!  hot  wall  by  the  action  of  the  nucleate  bubbles.  Gunther  (Cf.  Ref.  15)  performed 
siiiiilar  calculations  and  suggested  that  this  same  mechanism  of  induced  liquid  turbulence  may 
account  for  the  high  heat  transfer  rates  in  subcooled  liquids. 


B.  Bubble  induced  Liquid  Turbulence 

It  was  noted  in  the  discussion  of  bubble  dynamics  that  the  liquid  Su.Tooiidiag  the  bubble  wufi 
forced  away  from  the  heat-transfer  surface  during  the  growth  period  of  the  bubble.  For  vapor  bubbles 
forming  in  subcooied  watci,  the  liquid  was  Ors".™  back  toward  tit?  wall  during  the  collapse  period  of 
the  bubble  cycle.  The  resulting  liquid  flow  pattern  is  shown  in  Figure  89  (a).  The  same  liquid  flow 


rags  31 


5 

5 

t 

n 

£ 

j 

i 

i 


taut iWsfWi 


*ss 


1 


Memorandum  Ho.  20-88 


jaf  r ropij.1  k wn  L.QbvTGivjy£ 


puiieru  would  sc!  b“  fr-'mod  by  the  gea-vapor  bubbles  or  carbon  tetrachloride  bubbles.  Since  the 
latter  two  types  of  bubbles  detached  from  the  heat  transfer  surface  before  collapsing,  the  induced 
liquid  flow  would  be  as  shown  in  Figure  S3  (L),  U is  believed,  however;  'ha*  the  liquid  velocity 
induced  by  either  the  nondetaching  bubbles  nr  the  detaching  bubbles  may  be  related  to  the  bubble 
growth  and  collapse  velocities,  la  attempting  a correlation,  it  seems  reasonable  to  select  the 
average  bobble  velocity  (i„e.  max.  babble  diameter/bubble  lifetime)  as  the  typical  velocity  for  both 
the  degassea-waier  bubble  and  the  carbon  tetrachloride  vapor  bubble.  The  iocai  induced  liquid 
velocity  may  be  related  to  this  typical  velocity. 

In  addition  to  depending  on  the  bubble  velocity,  the  average  liquid  velocity  over  the  entire 
heat-transfer  surface  also  would  depend  on  the  distribution  of  the  bubbles  causing  the  liquid  velocity. 
In  order  to  include  the  effect  of  the  distribution  of  bubbles  on  the  induced  liquid  velocity,  a bubble 
population  parameter,  will  be  introduced,  so  that 

VL^  (Od>2^ 

• 

v^  « effective  velocity  cf  the  liquid  to  cool  the  wall 
£ - bubble  population  parameter 
f ■ bubble  velocity 

To  relate  the  effective  liquid  velocity  to  the  beat  flux,  we  may  try  the  uauai  dimensionless 
parameters  that  have  correlated  successfully  forced-convection  beat  transfer.  The  Nusselt  number 
is  related  to  the  Reynolds  and  Prandtl  number  as  follows: 

Nu  - CReaPrjS 

Substituting  the  effective  velocity  ana  bubble  radius  into  this  relation  results  in  the  following  ex- 
pression for  nucleate  boiling  heat  transfer: 


--C, 


/ PC*rm„Y/cl*Y  (r  T\ 

v- ri  U)  Vw'Tl) 


(16) 


where:  {£i)  has  replaced  the  free-stream  velocity  and  the  maximum  bubble  radius  has  replaced  the 
hydraulic  diameter. 

It  was  shown  in  the  introduction  that  the  principal  information  needed  for  the  design  cf 
heat-transfer  apparatus  operating  in  the  nucleate-boiling  region  are  the  peak-beat-flux  values  for 
a range  of  liquid  conditions.  As  a result,  it  would  be  important  to  rewrite  Equation  (16)  to  apply  at 
tbe  peak-beat-flux  point.  Tbe  Equation  may  be  rewritten  after  recalling  some  of  the  experimental ly 
determined  facts  about  nucleate  boiiing. 

It  was  seen  in  Section  IV -C  of  Puri  !,  that  the  heat  flux  in  the  nucleate-boiling  region  cotrid 
be  increased  appreciably  with  little  change  in  the  wall  temperature  or  bubble  velocity.  Thus,  if 
Equation  (It)  represents  the  beat  transfer  properly,  the  heat  transferred  by  nucleate  boiling  may  be 
increased  from  incipient  boiling  to  the  peak  heat  flux  with  the  major  change  being  the  population  of 
the  nucleate  bubbles.  The  population  parameter  £ which  relates  bubble  velocity  to  induced  liquid 
velocity  should  be  a function  of  (1)  the  ratio  of  average  distance  between  bubbles  to  ibe  average 
bubble  radius  and  (2)  tbe  Reynolds  number  based  on  bubble  velocity.  As  tbe  bubble  population 
increas^n.  a point  would  bo  reached  at  which  the  average  induced  liquid  velocity  would  be  a maximum. 
This  point  of  maximum  liquid  velocity  should  correspond  to  tbe  point  of  peak  ucat  flax.  As  s firs! 
annmyimatinn.  s!  seems  nermiKHihle  to  replace  the  popuiuiiou  parameter  by  the  product  of  the  Revnnld* 
number  to  soruss  power  and  a tactor  which  is  '■datively  insensitive  to  the  liquid  properties  aud  con- 
ditions at  the  peak  heat  flux  point,  Tbe  validity  of  this  assumption  may  be  verified  by  rewriting 
Equation  (16)  and  comparing  the  predicted  values  of  peak  heat  flux  with  the  experimental  data.  At  the 
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ooint  of  maximum  liquid  velocity,  it  oiems  permissible  to  replace  this  population  parameter,  with 
a factor  thai  should  be  approximately  constant  regardless  of  liquid  properties  or  coudiiioas,  i nerciore, 
iua  y wii tc  the  peak  k==t  flax  as 


(1?) 


where:  C2  has  replaced  the  product  Cj(£)aof  Equation  (16). 

The  three  constants  appearing  in  Equation  (17)  may  be  evaluated  by  choosing  the  values 
which  best  fit  the  experimental  data  for  distilled  water,  water-Aerosol  solution  and  carbon  tetrachloride. 
The  values  were:  C2  **  0.053,  a - 0.8,  fi  - i.0. 

The  calculated  peak -heat- flux  values  for  distilled  water  and  water=Aeroso!  solution  are  shown 
in  Figure  90.  It  is  seen  that  the  experimental  peak  "heat-flux  points  are  within  15  per  cent  of  the  esti- 
mated curve.  The  calculated  peak-heai-flux  values  for  degassed  carbon  ietra chloride  arc  plotted  in 
Figure  91.  The  experimental  data  again  are  seen  to  be  within  15  per  cent  of  the  estimated  curve. 

It  will  be  of  interest  to  consider  the  factors  responsible  for  the  trend  and  the  relative  positions 
of  the  peak-hcat-flux  curves  plotted  in  Figures  38  and  41  for  the  degassed  water,  Aerosol-water  solution, 
aerated  water,  degassed  carbon  tetrachloride  and  aerated  carbon  tetrachloride.  The  trend  and  relative 
positions  of  the  curves  are  discussed  in  the  foi'owing  Section. 


C.  Effect  of  Liquid  Conditions  on  Poak  Heat  Flux 

1.  Shape  of  Peak-Hsat-Flux  Curves,  As  the  hulk-liquid  temperature  wan  decreased,  the 
bubble  collate  velocity  increased  slightly  (Cf.  Sect.  IV-C-3),  and  the  physical  properties  of  the  liquids 
also  changed  slightly.  For  all  liquids  studied,  the  changes  of  bubble  velocity  and  physical  properties 
with  liquid  temperature  were  small  compared  to  the  change  in  the  factor  (Tj-  - Tj  ).  As  a result,  the 
peak  heat  flux  decreased  with  an  increase  of  liquid  temperature. 

2.  Relative  Positions  of  Pcak-Hcnt-Flux  Curves,  To  investigate  the  effect  of  surface  tension 
on  the  boiling  process,  one  per  cent  of  Aerosol  was  added  to  the  water.  The  photographic  study  showed 
that  the  bubble  velocity  decreased  to  about  A the  value  for  the  bubbles  growing  in  distilled  water.  The 
lower  bubble  velocity  reuulted  from  the  decrease  in  surface  tension  for  the  reasons  discussed  in  Section 
IV -F.  From  Equation  (17)  it  can  be  predicted  that  the  lower  bubble  velocity  should  produce  a lower  peak 
heat  flux.  It  in  seen  in  Figure  38  that  the  burnout  values  for  the  Aerosol-water  solution  are  consideralJy 
lower  than  for  the  distilled  water. 

The  carbon  tetrachloride  bubble  velocity  was  about  % the  distilled  water  bubble  velocity.  The 
lower  bubble  velocity  coupled  with  lower  thermal  conductivity  resulted  in  lower  peak-hest-flu*  values 
for  the  carbon  tetrachloride  (Cf.  Fig.  41). 

It  was  proposed  in  Section  1V-E  that  the  amount  of  gas  dissolved  in  a liquid  up  to  the  saturation 
condition  corresponding  to  the  wall  temperature,  would  be  expected  to  decrease  the  bubble  velocity  a 
snail  amo'j"*-  However,  increasing  the  concentration  of  gas  above  the  wail  temperature  saturation 
condition  would  be  expected  to  cevae  v significant  decrease  in  the  bubble  velocity  and  resuit  in  stagnant 
bubbles.  It  was  also  noted,  as  predicted,  that  the  difference  between  the  wall  temperature  and  the  satu- 
ration temperature  was  slightly  decreased  (Cf-  Fig.  39),  From  Equation  (17).  it  is  predicted  that  the  lower 
wail  temperature  and  smaller  bubble  velocity  would  result  in  a lower  peak  heat  flux.  It  is  seen  in  Figures 
38  and  41  that  both  aerated  water  and  aerated  carbon  tetrachloride  had  lower  peak-heat-flux  values  than 

- i*  i i i : . ; j„. 
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3.  Liquid  Pressure  and  Free-Stream  Velocity.  The  present  investigation  of  the  mechanism  of 
nucleate  boiling  heat  transfer  was  confined  to  a study  of  atmospheric  poo!  boiling.-  Althougn  no  data  were 
taken  to  study  the  effects  of  liquid  pressure  and  tree-stream  velocity,  itwiii  he  of  interest  to  consider, 
briefly,  the  possible  influence  of  these  two  factors  on  the  boiling  process. 
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!»  wai  seen  «n  Section  IV-C-3  that  the  bubble  velocity  would  be  expected  to  decrease  for 
an  increase  o f liquid  pressure.  As  a result,  the  heat-transier  coefficient  would  teed  ts  be  lower  »» 
high  pressures.  This  decrease,  however,  may  be  overcome  by  the  hitter  wall  temperature  and  re- 
sulting lower  liquid  film  viscosity.  Since  20  data  are  available  for  the  bubble  velocity  at  high 
pressures,  it  is  not  possible  to  check  the  validity  of  the  proposed  correlation  ia  Equation  for  a 
range  of  pressures. 

The  effect  of  forced  convection  02  boil*"®  heat  transfer  also  has  not  been  studied  in 
detail.  It  is  knows  that  the  nucleate  bubbles  slide  downstream  at  about  ^10  of  the  frec-stream 
velocity.  The  sliding  motion  may  be  visualized  an  coupling  with  the  bubble  velocity  to  produce 
additional  turbulence  in  the  liquid  boundary  layer.  Regardless  of  the  manner  in  which  the  free 
stream  velocity  and  the  bubble  velocity  combine,  it  would  be  expected  that:  (1)  an  increase  in 
thv  frrr  bucoui  velocity  would  increase  the  peak  heat  flux  values  (2)  for  free  stream  velocities 
much  lower  than  the  bubble  velocity,  the  peak  heat  flux  should  Le  insensitive  to  changes  in 
free-BUeaui  velocity  and  (3)  for  frec-stream  velocities  much  greater  than  the  bubble  velocity,  the 
peak  heat  flux  shonld  be  more  sensitive  to  changes  in  the  free-stream  velocity. 

These  general  predictions  have  been  verified  experimentally:  (1)  Figure  92  from  Reference 
15  shows  that  the  peak  h®ai  flux  increases  with  velocity  (2)  McAdams  found  for  low  free-stream 
velocity  that  the  peak  heat  flux  conelated  with  the  free-stream  velocity  to  the  l/3  power  for  ve- 
locities from  1 to  12  ft/sec  and  (3)  Gunther  found  for  high  free-Btream  velocities  that  the  peak 
heat  flux  correlated  with  free-stream  velocity  to  the  'A  power  for  velocities  from  5 to  40  ft/sec. 

D.  Conclusion 

The  agreement  between  the  calculated  peak  heat  flux  and  the  measured  burnout  points 
supports  the  assumption  that  the  principal  mechanism  of  heat  transfer  by  nucleate  boiling  is  not 
letent-heat  transport  but  rather  convective  heat  transfer  by  the  induced  liquid  velocity.  As  dis- 
cussed previously,  the  bubble  velocity  data  may  be  in  error  as  much  as  20  to  30  per  cent,.  However, 
the  predicted-heat-flux  curves  would  still  be  in  reasonable  agreement  with  the  experimental  data. 

In  view  of  the  small  heat  transfer  that  may  be  accounted  for  by  latent  heat  transport  and 
the  agreement  between  the  predicted  and  experimentally  determined  value"  of  peak  heat  flux  when 
using  Equation  (17),  it  seems  likely  that  the  induced  turbulence  in  the  liquid  boundary  layer  is  the 
basic  mechanism  by  which  heat  is  transferred  in  the  nucleate  boiling  region. 
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C,  Cj , C2 


• area  (cq  in.) 

■ velocity  of  sound  (in. /sec) 

■ specific  heat  (Btu/lb) 

“ cutinitoiia  (dimensionless) 

■ specific  heat  at  constant  pressure  (Btu/'lb-°F) 

• thermal  diffusivity  (sq  in. /sec) 

■>  outside  diameter  of  tube  (in.) 

« equivalent  diameter  (ia.) 

«»  voltage  drop  across  the  heating  tube  (volts) 

« acceleration  of  gravity  (in./sec^) 

■ mass  flow  rate  per  unit  area  (lb/sq  in. -sec) 

- total  heat  transfer  coefficient  (Btu/sq  in.-»ec-°F) 

- conduction  heat  tranafer  coefficient  when  radiation  is  negligible  (Btu/sq  in. -sec) 

■ conduction  heat  transfer  coefficient  when  radiation  is  appreciable  (Btu/sq  in.-sec-°F) 

- radiation  beat  transfer  coefficient  (Btu/aq  in.-sec-0  F) 

«*  current  (amps) 

* conversion  tactor(~5u/aec  'j 

\ watt  / 

» thermal  conductivity  (Btu/sq  lfi.-sec-°F/in.) 

- length  o{  heat  transfer  surface  (in.) 

- evaporation  rate  per  unit  area  (lb/sq  in.-sec) 

* Nuaselt  number 

» pressure  (lb/sq  in.) 

» pressure  of  gas  (psi) 

- pressure  inside  bubble  (psi) 

■ presntre  of  liquid  (psi) 

■*  pressure  cf  vapor  (psi) 

= Prandti  number  [Cp  uA  (dimensionless)] 

= heat  flux  (Btu/sec) 

« heat  transfer  by  conduction  (Btu/sec) 

«=  heat  transfer  by  radiation  (Btu/sec) 

- heat  flux  per  unit  area  transferred  to  the  test  fluid  (Btu/sq  in.-sec) 


/Btu/sec 


hr;:?  flux  per  uuit  area  transferred  to  the  stabilizing  fluid  (Btu/sq  in.-sec) 
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TABLE  I (Conn’d) 


^ ■ total  heat  flux  per  unit  area,^  — ^ ^ (Btu/sq  in.-sec) 


Q - 

r « 

t - 

ri  " 

n . 

r = 

Re  - 

t - 

T - 

Tf  - 

TL- 
T0  " 

T„.  - 

Tv  - 
Tw,  - 


T 

4 © ^ 


ATa  - 
AT  - 

ATS  - 

AT 
* sat 

AT  k * 

x sub 
u *■ 

V « 

VL  " 

W6  - 

x » 


y 

y. 

Vo 

• *• 

v+ 


heat  content  (Btu) 
radius 

bubble  velocity  (ft/sec) 
inside  radius  of  tube  (in.) 
electrical  resisiruce  (ohm) 
gas  constant  (ft/"R) 

Reynolds  number,  (dimensionless) 

tube  or  vapor  thickness  (jn.) 
temperature  (°R) 

arithmetic  mean  temperature  of  the  vapor  film  (°F) 
liquid  temperature  (°  F) 

maximum  temperature  inside  of  tube  wall  (°F) 
liquid  saturation  temperature  (°F) 
vapor  saturation  temperature  (°F) 
tube  wall  temperature  next  to  test  fluid  (°F) 

tube  wall  temperature  next  to  stabilizing  fluid  (°F) 
temperature  rise  of  test  fluid  in  annulus  (°F) 

TWa  ” TL  or  TWg~  "rL  ''**) 

temperature  vise  of  stabilizing  fluid  (°F) 

T , o'  t°tr\ 
iWA  “sat  ' * ' 

Tsat-T^P) 

vapor  velocity  in  direction  of  tube  axis  (in./sec) 
vapor  velocity  in  direction  of  tube  radius  (in. /sec) 
effective  liquid  velocity  (ft/sec) 
mass  of  gas  (!bu) 

direction  parallel  to  tub?  axis  (in.) 

* i » *.  /■■  / lb  \ 

vapor  turn  Doay  iorce  per  unu  massi  iu/ 


in.  /sec* 


direction  parallel  to  tube  radius  (in.) 
distance  from  inside  tube  wall  to  point  of  maximum  temperature  (in.) 
distance  from  outside  tube  wall  to  point  of  maximum  temperature  (in.) 
dimensionless  distance 
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TABLE  I (Conl’d) 


t *•  accomodation  coefficient 
( y m emissivity  of  tube  wall  (dimensionless) 

,,  „ emmisivitv  of  liquid  (dimensionless) 

Is 

A - latent  heat  (Btu/lb) 

■ liquid  density  (!b/cu  is.) 
n - vaDor  density  (lb/cu  in.) 

iv  • 

a _ Stefan  Raltzmann  constant  (Btu/sq  in. -sec*  R) 


T0  “ 

e - 

e - 

v-  * 
V2  - 


.«•  -I 


surface  tenaiou 

shear  stress  at  wall  (Ib/in.  ' 
time  (sec) 

sk+  2x 

Ji.. 

VA  t// 

9 \ 

viscosity  (lb-sec/ia.  ) 

<£.+  <? 

^v2 

V — - / 

bubble  population  parameter 
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TASLc,  iii 

EXPERIMENTAL  DATA  AND  CALCULATED  VALUES 
USING  NITROGEN  GAS  STABILIZER 


Press.  = 16  peia 
Subcooling  - SCPF 


Vel.  ■ 1.1  ft/sec 


No. 

T 

1 

Amps 

E 

Volts 

(1) 

\a/T 

Btu/sq  in.-sec 

AT  A 
°F 

Btu/sq  in.-sec 

n 

v a/g 

Btu/sq  in.-sec 

twa 

O VO 

r 

10 

150 

3.1 

0.19 

0.9 

0.04 

u,.  15 

X 7U 

11 

1S6 

3.8 

i»3 

0.06 

0,23 

210 

12 

228 

4.7 

0.44 

1.9 

G.09 

0.35 

231 

13 

r\  . 9" 

5.2 

0.52 

2.5 

0.12 

0.40 

246 

14 

286 

6.1 

0.72 

5.1 

0.25 

0.47 

252 

15 

360 

7.6 

1.11 

11.7 

0.57 

0.54 

265 

16 

436 

9.2 

1.63 

20.2 

0.99 

0.64 

285 

17 

490 

10.8 

2.15 

21.5 

1.05 

1.10 

376 

18 

500 

11.2 

2.27 

19.5 

0.95 

1.32 

420 

19 

540 

12.9 

2.82 

13.5 

0.66 

2.16 

545 

20 

544 

13.1 

2,88 

8.2 

0.40 

2.48 

615 

Ol 

4<A, 

562 

14.0 

3.19 

5.2 

0.26 

2.93 

690 

92 

640 

17.6 

4.55 

2.2 

0.11 

3.44 

970 

23 

775 

22.7 

7.13 

2.5 

0.12 

« > 
i lUi 

1 POA 

1U4V 

24 

786 

23.0 

7.34 

2.9 

0.14 

7.20 

1560 

25 

812 

24.5 

8.05 

3.0 

0.15 

7.90 

1710 

Press.  » 16  psia 
Subcooling  =.  50°F 

Vel.  - 5 ft/sec 


26  1 

162 

3.4 

0.22 

1.9 

0.09 

1 

0.13 

190 

27  ' 

ZiZ 

5.7 

A 

V#l*J 

5.6 

0.27 

0.36 

235 

28 

332 

7.0 

0.94 

9.8 

0.48 

0.46 

265 

29 

456 

9.6 

1.77 

26.7 

1.30 

0.47 

Z iu 

30 

530 

11.6 

2.50 

1 1 o O 

OU«7 

1.90 

0.60 

ooc 

31 

540 

11-9 

'2.59 

40.0 

1.95 

0.64 

345 

32 

545 

13.0 

o a*r 

<--  • O i 

21.5 

1.05 

1.82 

| 537 

33 

550 

13.2 

2.94 

15.4 

0.95 

1.99 

540 

34 

555 

13.9 

&12 

12.3 

0.60 

2.52 

630 

585 

14.6 

3.46 

6.5 

0.32 

3.14 

765 

36 

6G0 

15.6 

3.79 

3.1 

0.15 

3.64 

«l  A 

OiU 

37 

685 

19.8 

5.50 

2.0 

0.10 

5.40 

1350 

38 

700 

21 .0 

5.95 

2.2 

0.11 

5.84 

1450 

39 

741 

l 

22.2 

L 

6.65 

2.2 

I 

°-U  J 

6.54 

ioGO 

rage  3? 


: n 

I 

i 

i 

} 


- i • • 


-rVi,,..  ' - 


TABLE  III  (Cant’d) 


Press,  - 16  psia 
Subcooling  ■ lOC^F 

Vs!,  ••  1. 1 ft/gAr 


Preset*  - 60  psia 
Subcoolii'g  ~ 60°  F 

Vei.  ■ lil  St/ sec 


17.4  | 0.85 

038 
030 
0.13 
0.12 
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APPENDIX 

ANALYSIS  OF  FILM  BOILING  ON  A VERTICAL  HEATING  SURF  ACE 


A rigorous  analysis  of  the  fiim-boiiing  heater  anaier  nrct'linuism  is  beyond  the  scope  of  this 
experunenta!  investigation.  However,  a simple  analysis  may  be  carried  out  and  ihe  results  com- 
pared with  the  experimental  data  by  employing  the  following  assumptions: 

1.  The  vapor  flow  is  steady  with  physical  properties  independent  of  temperature. 

2.  Vapor  boay  forces  are  negligible  compared  with  viscous  pressure  and  inertia 
forces. 

S.  1 he  change  in  vapor  velocity  along  the  tube  is  small  compared  with  the  radial- 

velocity  gradient  (i-e,  du/dx  « du/dy).  It  is  also  assumed  ihai  u(du/dx)  « \{da/dy). 

4-  The  evaporation  rate  is  constant  along  the  tube  (i.e,,  m - p'.'r  m(x,  y). 

5.  The  vapor  density  is  independent  of  the  radial  (y)  position, 

[i.e.,  pv  *-  pv  (x)  py  (y)]. 

6.  The  pressure  drop  along  the  tube  du-  to  liquid  flow  is  much  lens  than  the 
buoyant  forces  (i„e»,  dp/dx 

7.  The  water  velocity  at  the  water  vapor  interface  is  small  compared  with  the 
average  velocity  of  the  vapor  in  the  film  (i.e.,  the  vapor-film  velocity  is 
substfintieliy  zero  at  the  tube  surface  and  the  water  interface). 

8.  Heat  it*  transfened  across  the  vapor  film  from  the  wall  to  the  water  by 
conduction  and  radiation  only. 

9.  All  the  heat  leaving  the  tube  is  used  to  vaporize  the  liquid.  (This  assumption 
neglects  the  heat  that  may  be  transferred  to  the  bulk  of  the  liquid  as  well  as 
the  heat  that  may  be  absorbed  by  the  vapor  to  increase  itn  temperature  above 
the  saturation  point.) 

10.  The  vapor  film  is  smooth. 

11.  The  vapor-liquid  interface  is  at  the  saturation  temperniuio.  The  genera! 
momentum  equation  in  the  x direction  may  be  y/ritten  as 


r ““  "r  <r  v 1/  ^ r"7 9 / » i 

+ — X + 73  — + uv*u 

g dd  dx g dx 

Lisins  assumption  1 through  5 and  neglecting  the  effect  of  evaporation  on  tl*  momentum 
change,  this  equation  may  be  simplified  for  the  vapor  film  as  follows: 

d2*’,  1 dp 

- 0 (A -2 

A. 2 3 dx 


By  assumption  6 


• z 

Q u PL 


+ ^ - 0 


Integration  of  liquation  (A-3)  gives 


u •=  Cj  Cj;  ~~~r 
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The  boundary  conditions  given  by  assumption  7 are 

u « 0 aty  - t and  y - 0 

Then 


. r ..  ..2 

n.  / _ /_ 

p L 2 2 


(A-5) 


The  velocity  may  be  eliminated  and  an  expression  obtained  for  the  film  thickness  t by  the  continuity 
requirement 

x t 

r r 

! .i-TVK 


j mi?  Ddx  **  J pyn  Dudy 
0 0 


(A- 6) 


where  x - 0 at  the  start  of  the  heating  section. 

For  m and  p„  iiid-pendent  of  x and  y,  respectively, 


/ udy 
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The  velocity  may  new  be  eliminated  by  substituting  Equation  (A-5)  in  Equation  fA-7)  and  integrating 


mx  p^  td 
o..  u 12 


(A-S) 


For  the  case  where  the  radiant  heat  transfer  is  negligible,  assumptions  9,  10  and  11  may  he  used  to 
eliminate  the  evaporation  rate  fron  Equatiou  (A-8) 


3-k 


at. 


sal 


- mA 


or 


kllsii 

tA 


Then  Equation  (A-9)  and  (A-8)  lead  to 


t*  - 12— — 1 — 1 ±- 
A PLPy 


(A-9) 


(A-10) 


or  defining  a heat  transfer  coefficient  for  the  heat  conducted  across  the  vapor  film  at  point  x 


Then 
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Equation  (A-ll)  represents  the  heat-transfer  cocliicient  at  the  p 
average  heat-transfer  coefficient  for  a heating  length  L 


usiiioi:  v.  •"Ktair,  ths 


- r ^ 

CX  dL 


/°T  , -ix 

j ApLpv  k3  dx 
/ [l2ATBatxp]  L 


r „ 1% 

*P lP yt*  ! 


KL“V“,.-  (A-12) 

L^“  1 sat  ^Pj 

where  A Taat  is  the  average  value  between  x « 0 and  x « L. 

It  is  not  surprising  to  note  that  Equation  (A-12)  with  the  exception  of  the  factor  12  is  the 
same  as  the  heat-transfer  coefficient  for  film  condensation.  The  constant  factor  12  diffei-*  because  of 
the  difference  in  boundary  condition  at  y • t (Cf.  Ref.  6 and  7). 

The  heat  transfer  by  muiaiiou  may  be  expressed  by  the  conventional  relation 


R AT 


!>WA+460)4-(T8at  + 460>* 


1 t i 
a at  — + — ±-  ~ 1 I 

e*  e, 


(A-13) 


This  relation  is  valid  if  all  the  heat  is  absorbed  in  the  water  near  the  vapor-liquid  interface. 

Since  the  radiant  heat  transfer  affects  the  thickness  of  the  vapor  film,  the  conduction  and 
radietion  are  inextricably  related.  An  appreciable  deviation  from  a rigorous  solution  is  neccaeary  to 
account  for  the  radiation  in  a simple  manner. 

Equation  (A-9)  becomes 


& _ k A Taat  ^ ^ 3R 
tA  dr 


and  Equation  (A-10)  is  modified  to 


i* » 12  qc  + 


?vPLA 


The  heat-transfer  coefficient  tit  point  x ie 


Vi  qc  \*MPvP  V 


t \qc-qn/  \l2ATsatx 


or  the  average  heat-tranafer  coefficient  over  a distance  L is 


hc  - % ( 
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For  the  range  of  variables  encountered  in  this  study,  the  factors  (q^/q^.  + qR)^  varied  less 
than  20  per  cent  from  unity.  The  error  in  totui  neai  transfer  that  would  result  from  neglecting  this 
variation  would  be  less  than  10  per  cent.  The  average  total  heat-transfer  coefficient  for  the  length 
I . may  be  written  as 


hc  + hR 


73 
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A comparison  of  the  experimental  data  and  the  heat-transfer  coeffic'cnt  calculated  by 
Equation  (A-14)  is  shown  in  Figure  32  where  the  calculated  heat  flux  is  slightly  lower  than  the 
experimental  data. 

A urtuilcu  study  of  the  vttpur  film  was  beyond  the  scope  of  the  preseni  investigation. 
However,  5?  is  interesting  to  compare  the  film  thickness  calculated  from  Equation  (A— 10)  with  the 
film  thickness  estimated  from  the  photographs.  By  projecting  the  16-mm  picture  so  that  the  image 
was  enlarged  100  times,  it  was  possible  to  estimate  the  film  thickness  using  the  0.250-inch-iiameter 
heating  tube  as  the  base  scale. 

The  vapor-film  thickness  was  calculated  from  Equation  (A-10)  and  the  data  for  test  number  39 


i » 


12 


kA!«aL^ 

Pi 


n & 


.■jd 
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'1 

i 


12  r,.16xl0^^80)3(1.81xl0-9)- 
L 1000(2.29  x 10-5)  (3.62  xl0“2)  . 


- 0.0136  in. 


Th#.  vapor  film  in  the  photograph  for  test  39  was  approximately  0.020  inch  thick.  The  32  per 
cent  diffetvnce  between  the  calculated  and  measured  values  is  less  than  the  probable  error  in  esti- 
mating the  vapor-film  thickness  shown  in  the  photographs. 

An  interesting  point  in  the  region  of  film  boiling  heat  transfer  is  the  condition  at  which  the 
smooth  vapor  film  transposed  into  a rough  film.  !t  is  reasonable  to  suppose  that  this  transition 
would  occur  when  the  vapor  in  the  film  changed  from  laminar  to  turbulent  flow.  Using  the  data  from 
test  23  and  an  arithmetic  mean  temperature,  the  Reynolds  number  based  on  equivalent  pipe  diameter 
was  220.  The  vapor  film  for  this  test  was  rough.  The  critical  Reynolds  number  for  turbulent  flow  in 
a pipe  with  forced  convection  is  approximately  3000.  In  view  of  this  high  critical  value,  the  Reynolds 
number  of  220  obtained  for  ibis  test  should  have  been  too  small  to  produce  turbulent  flow. 

The  Reynolds  number  based  or.  distance  from  the  leading  edge  for  flow  over  a flat  plate  was 
approximately  14,000  compared  with  the  conventional  critical  value  of  80,000.  Here  again  i:;=  value  of 
Reynolds  number  is  well  beiow  the  conventional  critical  value. 

Correcting  both  of  these  Reynolds-numbur  criteria  for  non-isovhermal  flow  with  mass  addition 
probably  would  not  account  for  the  large  discrepancies.  It  is  possible  that  the  rough  vapor  film  resulted 
trom  orncr  exiero&i  iuclCm  o o iiv- u v3r»«t  10 »is  in  water  flow  or  pressure  surges  from  the  pump. 

Gne  more  rough  quantitative  comparison  may  be  made  between  the  experimental  data  and  the 
calculated  values.  The  mean  vapor  velocity  for  steady,  uniform,  viscous  flow  with  cccrlcretj-g 


\T  dp 
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Substituting  the  relation  for  t from  Equation  (A-10)  gives 
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/kATsa-  *f>l. 

UmV  12  Apvp 


/ 8.14  x 10-7  (1200)  3 (. 
y 12  (970)  (1.21  x 1(T5)(3.’ 


356  in. /sec  = 29.6  ft/sec 


.0361) 

78  x 10“9) 


Tlie  estimated  mean  velocity  obtained  from  measuring  the  photographic  data  was  roughly  25  ft/sec. 


it  in  concludes  last,  over  mr.ge  oi  ver,“!>!?s  investigated,  the  derived  equation  wi!!  give 
satisfactory  engineering  estimates  of  the  heat  transfer  in  the  complete-film-boiling  region  for  a vertical 
heat-transfer  surface.  However,  a more  rigorous  analysis  should  be  made  and  experimentally  verified  to 
give  better  understanding  of  the  process.  Further  study  of  the  complel.e-film-boiling  region  should 
include  measurements  of  the  vapor-film  thickness,  velocity,  and  stability  characteristics.  With  slight 
modifications,  the  equipment  tie  scribed  in  this  report  would  be  suitable  for  this  detailed  study  of  the 
complete-film-boiling  region. 
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Figure  38.  Peak-Heat-Flux  Values  for  water 
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Bubbie  Radius  vs  Time 
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Bubble  Radius  vs  l ime 
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Figure  53 

Bubble  Radius  vs  Time 
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Bubble  Radius  vs  Time 
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Figure  57.  Hubble  Radius  vs  Time 
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Bubble  Radius  vs  Time 
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Bubble  Radius  vs  l ime 
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Bubble  Radius  vs  Time 
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Sketch  of  Nucleate  Bubble  in  Degassed  Water 


Figure  79.  Temperature  Profile  in  Bubble 
Film  During  Growth  Period 


Figure  80.  Temperature  Piofile  in  Bubble 
Film  at  Maximum  Radius 
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Figure  81.  Excess  Wall  Temperature  vs  Liquid  Pressure 
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Figure  86b.  Sketch  of  Nucleate  Bubble  in  Degassed  Carbon  Tetrachloride 
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